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ABSTRACT 


The temperature dependence of the longitu- 
dinal sound velocity in single crystals of hcp He4 
was measured from .75 K to the melting point at 86 
bar and 120 bar with a 5 MHz longitudinally-cut 
quartz transducer. In the high temperature region 
the velocity obeys a simple power law. Itsdirection 
of change is contrary to that predicted by a reduced 
equation of state. At ee. Lona = 6/20 a sharp 
anomaly or "knee" is observed. The resulting plateau 
(the velocity changes only slightly in this region) 
extends to lowest temperatures. We feel this plateau 


is evidence of the long sought coupling between first 


and second sound velocities. 
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CHAPTER 1 


INTRODUCTION 


Solid helium was first studied in the hope of 
finding the simplest and most ideal of solids. These 
early experiments indicated that, in gross respects, 
helium was a solid like all other solids. However, 
detailed quantitative studies of the solid produced 
results remarkably different from the theories of the 
day. 

The helium atom is quite light, and in the weak 
van-der-Waals potential well of its neighbors in the 
solid, has a large zero-point motion. This quantum-like 
behavior has led to investigations of helium and similar 
materials under the general name quantum solids. All 
modern theories of solid helium directly take into 
account the quantum nature of the solid and have explained 
many observed “anomalies", although many unanswered 
questions remain. 

There are two stable isotopes of helium, He? 
and He’. In the liquid state the two isotopes behave 
drastically different. This is due to the fact that 
He4 has even spin (2 neutrons plus 2 protons plus 2 
electrons) and thus obeys bose-einstein statistics, 
whereas He? has odd spin (1 neutron plus 2 protons plus 
2 electrons) and thus obeys fermi-dirac statistics. 


These statistical effects come from overlap of the 
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individual wave functions. In the solid, the particles 
are highly localized with very little overlap. The over- 
lap, however, is large enough to lead to some quantum 
effects such as spin exchange in He? and short range 
correlation in both isotopes. 

All our experiments were performed on He* 
because chemical and isotopic purity can be obtained 
using simple and inexpensive techniques. Transport 
quantities, such as second sound and thermal conductivity 
are influenced dramatically by small concentrations of 
impurities. In solid helium, second sound can be 
propagated in a lightly damped mode while the thermal 
conductivity exhibits a large umpklapp maximum and a 
well defined phonon-Poiseuille flow region. This is 
experimental evidence that samples of He# low in 
impurities can be produced. 

The experiments reported on here mark the 
second stage of work to determine the temperature, 
pressure, and frequency dependence of the attenuation 
of first sound in #e*, Although some qualitative 
observations were made, an accurate determination of 
the attenuation proved too difficult with the present 
state of technology. 

However, measurements of the temperature 
dependence of the velocity showed anomalous behavior. 
This anomalous behavior appears to be evidence for 


coupling between first and second sound. The main 
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purpose of this work is to report on these findings. 
This anamalous bchavior has been seen before in liquid 
ne4 (Whitney and Chase (1967)) and predicted as a 
coupling between first and second sound by many authors, 
including Niklasson (1968, 1970, 1971) and Gurevich and 
Efros (1966, 1967). 
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CHAPTER 2 


THEORY 


So much work has been done on the theory of 
solid helium on both lattice dynamics and sound 
velocity that I will mention here only some of the 
review articles recently published. 

Classical lattice dynamics are treated exten- 
Sively in several books. I have relied mainly on 
the paper by Musgrave (1954) and the series Physical 
Acoustics edited by Mason (1965). Lattice dynamics 
of solid helium are well covered in the review papers 
by Werthamer (1969) and Guyer (1969). The latest 
experimental data will be found in Wilk's book (1967), 
Keller's book (1969) and Trickey et al (1972). The 
temperature dependence theories were obtained with 
consultation to Bhatia's book (1967) and the series 


by Niklasson (1968, 1969, 1970, 1972). 


2.1 Classical Lattice Dynamics 


We shall present here a short introduction to 
classical lattice dynamics based on Ziman (1965). 
Although this formulation is not sufficient to solve 
the special problems of solid helium, the initial steps 


and broad treatment are identical. 
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If we have a lattice with identical atoms of 
mass m at equilibrium lattice sites ~% and where u 
is the displacement of each atom from that site, then 


the kinetic energy of the crystal is T 


u : (2.1-1) 


We will now write potential energy in a Taylor 


. . . = 
series expansion in terms of the u_ so that 


(25 b=2) 


where j represents the cartesian components of u. If 
we truncate the series after three terms (as above), we 
call this the harmonic approximation. Note that the 
first term is unimportant, and the second term vanishes 
for small oscillations about the equilibrium position 
in the bottom of the potential well. 

If we use Lagrange’ formulas we can write the 


equations of motion as 


2 + 
i =~) j ar Ms 7 ] = : (2.1-3) 
£'5' on oe 

ae | 
If now we let oi represent the elements of the 


dynamical matrix G, then 
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Assuming translational invariance, we must have 


(201-4) 


- Tete 
GJJ = GJ (h) where i = L' - L (2.1-5) 
gh 
hence 
: ae oat 
One ee, te. ee ome (2.1-6) 
£ oe +h 
hj 


These several equations must satisfy the Bloch condi- 


tion 


u_ tt) = etd- u(t) = @id- u, (t) (2.1-7) 


hence, with the origin as arbitrary 
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We now assume a plane wave solution of the form 
ser ge (2.1-9) 
and we get in (2.1-8) 
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The above equation has a solution if and only if 
the secular equation (the equation of the determinant of 
the set of equations' coefficients set to zero) has a 
solution. This secular equation has 3n solutions where 
the n represents the number of particles in a unit cell 
and the three represents the three modes or polarizations 
of vibration (two transverse and one longitudinal). The 


sound velocity cJ (q) is given by 
ie) ee os hae Cr (2.1-11) 


We usually only use the term 'sound' in the limit 


w + O. In this region 


where j represents the mode. The dispersion in snitct’ 
is negligible until one reaches frequencies of order 
10? Hz. 

We note, also, that the elements of the dynamical 
matrix are closely related to the elastic constants in 


the long wave length limit (w + 0). This will be further 


dealt with in section 2.3. 


2.2 Lattice Dynamics of Solid He4 
The classical theory of lattice dynamics fails 


dramatically for solid helium. 
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If one chooses the interaction potential 
between two helium atoms to be a Lennard-Jones of 


the type 
Aes pert oti ened CP eo hel (2.2-1) 


and derives e and o from the virial coefficients of 
He® gas, the classical theory of lattice dynamics 
predicts a molar volume of 10 ef? (too small by a 
factor 2), a compressibility too small by a factor 30, 
and a sound velocity too large by a factor 4 when 
compared to hcp He at 0K and 26 bar the minimum pres- 
sure of solidification! Even worse, when De Wette and 
Nijkoer (1965) attempted to calculate the phonon fre- 
quencies and sound velocities using the aforementioned 
theory, they calculated imaginary frequencies for molar 
volumes greater than 12 cm? /mole (even though the 
lattice is stable up to molar volumes of 21 cm? /mole) ! 

On the other hand, thermal conductivity, elas- 
ticity and heat capacity measurements showed helium 
to behave quite normally. 

The major cause of these deviations from classi- 
cal behavior is the large zero-point motion of solid 
helium atoms. The attractive potential between helium 
atoms is the relatively weak Van der Waals interaction 
which produces a very narrow potential well. Estimat- 


ing the zero-point motion from the indeterminacy 
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principle yields for its kinetic energy 


at RYN Fee al 


z om 3 (2.2-2) 


where a is the lattice spacing. 

This large zero-point motion, combined with a 
narrow potential well, means that the mean position 
of the helium atom is not at the bottom of the well, 
but displaced to increase the lattice spacing. This 
means that the molar volume will be larger than 
expected as will the compressibility. 

Figure 1 shows the position of a helium 4 atom 
in the well of its nearest neighbors. The cross- 
hatched area represents an approximate wave function 
for He’, Notice that the He4 atom sits at a relative 
maximum in the well which accounts for the imaginary 
frequencies obtained by any classical calculation. 
Thus, the wave function has two peaks corresponding 
to the minima of the well. One can think of the atom 
as moving in a spherical shell about its mean position. 

The motion of any atom about its site will 
change the position of the potential well of its 
neighbor. Thus the neighbors motion will be correlated 
in such a way that the wave functions overlap as little 
as possible. We call correlations of this kind short 


range correlations. 
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FIGURE 1 He# atom in potential well of nearest 
neighbors (Guyer (1969)). 





ee eae wae | 
(IIIA) OV, 
a ° a 7 


at 
7 Of 





12 


If phonons exist in a solid, they are evidence 
of long range correlations in the solid. Any theory 
must therefore take into account the large zero-point 
motion, short range correlations and long range corre- 
lations if it expects to explain the behavior of solid 
helium to any accuracy. 

There are several theoretical approaches to the 
problem of lattice dynamics in helium. In general 
they can be divided into theories that use a single 
particle approach and those which use a many-body 
approach. 

The single particle picture (Nosanow (1964), 
Fredkin and Werthamer (1965)) assumes a Jastrow type 


wave function (Jastrow (1955) of the type 


y(ry---4)) = : o,(r, -S) f(r, - zr 


k? (2, 2-3) 
where >, are the single particle wave functions loca- 
lized at lattice sites a while £(r- as describes 
the correlation between particle 3 and K. We know 
something about the form of the correlation function 
f(r). For example, at large r it must be unity while 
at r=0 it must vanish because of the hard core. 

Also there should be a maximum in f about where the 


minimum in the Lennard-Jones potential V(r) occurs. 


A simple analytic form for f(r) is thus 
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HCP He4 
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Internal energy as a function of molar volume 
comparing theory and experiment (single par- 
ticle Nosanow (1966) and Hetherington et al 
(1967); collective - Gillis et al (1968). 
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f(r) = e CV") (2.2-4) 
where C is a variational parameter. Although there 
are many possible forms for f(r) consistent with the 
above criteria, they must be easily subject to varia- 
tional analysis or they are too difficult to use. 
Assuming an initial value for C, we can cal- 
culate the ground state energy Eo: If h, is the one 
particle hamiltonian and h, the two particle hamil- 


tonian we can write Eo in a cluster expansion 


<j{h, (5) | 53> 1 ) sheet ea <j{h, (5) |3> 
jk 


j <j|j> <jk|jk> <j|i> 


eee ly (Seu . 
<k|k> ane <jkz|jke> 
- two body terms + one body terms| (2.2-5) 


where j is the ¢hh 


particle wave function of jy. 

We truncate E after one term, minimize the equa- 
tion for E to get E, and then solve for |j> =. 

It turns out that $, is roughly a spherical 


gaussian (as for the ground state of the simple harmo- 


nic oscillator) : 
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a AS 


os =e (2.2-6) 
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where one minimizes E. for a given C. One then self- 
consistently minimizes C and A to obtain the lowest 
Eo: 

,.One can then take the cluster expansion to 
higher terms in E, although to my knowledge the 
expansion has not been carried past three terms. 

Using many body techniques, Fredkin and 
Werthamer (1965) have shown how to calculate the 
phonon frequencies from these single particle wave 


functions. 


SR ene ty eae ta Ga re <0|Vog¢(x-x' th) | 0>0 (4) 


at T = 0, where Vorfe is 
TEL: ee ee V-vi(r) Je ckY fF) (2.2-7) 
ote " 


and V(r) is the Lennard-Jones potential. This above 
calculation is the same as using VoFe in a Hartree 
calculation and taking no correlations into account. 


Many calculations have been done, but so far no cal- 
4 


culations have been done for hcp He 
In the theories emphasizing collective phenomena, 

one usually examines the response of a crystal to an 

externally applied disturbance. These so-called self- 


consistent theories (Hooton (1958), Bernades (1960), 
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Boccara and Sarma (1965), Raninger (1965), Koehler 
(1966, 1967, 1968), Choquard (1967), Horner (1971)) 
consider the phonons as the basic co-ordinates of the 
solid. Essentially the result of these theories is 
that the spring constant between a pair of atoms is 
the second derivative of the interaction potential 
averaged over the motion of the atoms. Short range 
correlations must be added in later. These theories 
appear much better suited for numerical analysis. 

If we say that the hamiltonian of our solid can 


be approximated by 


Ae BT + ST Sin: (2.2-8) 
: bok ) ag 
i afi 


and further approximate ‘# by Hy 


Action, . 
y i335" y3 


ype 3 eit 


Ni] FR 


2 
* fi 2 
be ET ak aa 


where the G's are the spring constants and will be 
solved for variationally. The energy is minimized 


(Koehler (1968)) by 


Gaia <o|v2 v2 |O> (2.2-9) 
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and the ground state eigenfunction is given by 
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Hence instead of averaging with respect to the single 
particle wave functions, we here are averaging with 
respect to the phonon ground state. 

A third approach (Guyer (1969)) treats both 
short and long range interactions equally by using 
a Rayleigh-Schroedinger perturbation approach. This 
theory offers a fair description of He? in lowest order 
with a well defined method for going to higher orders. 


Take the hamiltonian of the system to be 


ig oe ol 8 t Y 
where vi(ij) = v(ij) - Pr sad) and v(ij) is the pair 
interaction and Base)? is the interaction between this 
pair and the rest of the particles. One is allowing 
a set of one, two, three or more particles to interact 
with each other in the Hartree field of the other par- 


ticles 


Ao o(1---N) = EL 6, (1---N) (2.2-12) 


where @, is the ground state wave function. 
The Rayleigh-Schroedinger perturbation (Brueckner 


(1950)) expansion for the energy shift is 
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which is identical to solving the equation 


Pilla V(ij) 1955 (1---N) = £55055 (1---N) : (2.2-15) 


This is just our hamiltonian including two body terms 


only. To include three body terms 
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and subtraction occurs because these three terms al- 
ready appear in the two body term. 

In essence, the one body term is the self- 
consistent phonon theory with no short range correla- 
tions. The two body term is the self-consistent phonon 
theory with two body correlations. This theory then 
is able to treat long and short range correlations in 


a coherent systematic manner. 


2.3 Wave Propagation in Ideal Elastic Solids 


The elastic constants c. describe the stresses 


ZIKL 

O55 in Vawious drecefons/inythe crystal in terms of 
strains Erg applied in various other directions 

ee i L Ci ike" ok ; (2.3-1) 

in es 

55 is the force acting on plane i in direction j and 
the strain Ey is defined in terms of the displacement 
> 
u as 

e,, MS a ae (218-2) 
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Equation (2.3-1) is simply Hook's law. 
The equation of motion for an elementary cube 


is thus 
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Phonon dispersion curves in the Basal plane 
for the collective picture as calculated by 
Gillis et al (1968) compared to experimental 
data of Lipshultz:et al (1967). 
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This has a plane wave solution of form 
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We now may determine the sound velocities by solving 
the secular equation for q. 

If we now apply symmetry and thermodynamic krqu- 
ments for hexagonal crystals, eTakL has only 5 independent 
components Cry Shoe C31 C337 C,,(where we have used 


the shortened notation of Voigt - =c 


evi eee 1122 12° 
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Several authors (Zener (1936), Gold (1950), 
Musgrave (1954)) have explicitly solved the secular 
equation for hexagonal symmetry. They find velocity 
surfaces (which are completely determined by the elastic 
constants and the density) have cylindrical symmetry 
about the "c" axis and hence the velocities can be 


expressed in terms of one angle "y" the angle between 


propagation and the "c" axis. They find: 
% 
V, = (cy/p) 
% 
Va Ny (7) 
c 29-8) 
Wik (c, 7/0) (2. 
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FIGURE 4 


c-axis 





Velocity as a function of orientation accord- 
ing to a collective picture cajculation by 
Gillis et al (1968) for hcp He* at 80 bar. 
The curves have cylindrical symmetry about 
the 'C'-axis. 
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For the special case of propagation along the "c" 


axis (y = 0) 
Tmat faa 
(2.3-9) 
Se tn oe AA 
and in the basal plane (y = 90°) 
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Velocity surfaces are plotted in figure 4 accord- 
ing to Gilles et al (1968) for hep He4 at 80 bar. 

Now, in general, sound velocities do not propagate 
perpendicular to the transducer, but may, in fact, deviate 
quite dramatically from a normal direction. Figure 5 shows 
the calculated deviation for the elastic constants cal- 
culation by Gillis et al (1968). The method of determining 
these deviations is given in an appendix in the form of a 
computer program. 

For purely transverse waves, the disturbance must 
be perpendicular to the direction of propagation. For purely 
longitudinal waves, the disturbance must be along the direc 
tion of motion. For hexagonal crystals, only one of the 
transverse modes is a pure mode, the other two being mixed. 
That is,the disturbance moves at an angle to the direction 
of propagation that is a function of orientation. The two 
modes are always perpendicular to one another at any 
Orientation and one is predominantly longitudinal while the 
other is predominantly transverse. We call these modes 


quasi-longitudinal and quasi-transverse respectively. 


2.4 Related Quantities 
If the elastic constants of a crystal are known, 


then several other quantities can be determined 


Debye Temperature 85 
The Debye temperature at 0k is (Alers (1965)) 
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FIGURE 5 Propagation direction schematic. 


26 


as 





ye>ubenn1t 





erpliwe lotey12 
ire 


_ 7 : oa 


27 


where 


h = Planck constant 

k = Boltzmann constant 
N = Avogadros number 
/ = Molar volume 


A computer program in the appendix for both 
cubic and hexagonal crystals is given. The On cal= 
culated in this way should be identical to the 
calorimetric Debye temperature at temperatures below 


0,/50. 


Compressibility 
For hexagonal crystals, the isothermal compres- 


sibility is (Nye (1964)) 


Ss Ss Ss Ss Ss 
C773 gy tolyoSagi7 2¢73 


One must remember, however, that the compres- 


sibility calculated from elastic constants is adiabatic. 















\ 
é vb Ob ‘f 7 o~. ‘ t=) 7 . 
T > 23 i3 3 a 
toodenod AohaelI = dd’ 
tosjenoo ansmstiog = a 
zedmyvn zorbspévA = 4- 
amuilov yeloM = PAS a 
‘Joa 101 xifmegqs sis af msitpoxig tetuqmoo A “i 
“+59 0 ot? .nevip ef afegeyim [snopsxed bak sidgo 
ait 03 lsottnebi sd Bluote yew atd> ai bedaelgs 


* 


= 


J ? 
wOleG ao Tus ategmDs Js stude1aqmed sydad cixzemtzolse 


-02\,0 = 


i. 


yilidtasergmod 


“sesi1gmo> laemisitoet odd ,alstjeyio Isnopsxed xoO¥ 


_((aGl) oy) 2b yselidde 





28 


The difference between these two values is discussed 


thoroughly in sections 2.5 and 2.6, but is of order 


where Cp and cy are the heat capacities at constant 


pressure and volume. 


Phonon Spectrum 


9 


In the long wave length limit (w << 10° hz), the 


sound velocity is equal to the phase velocity 


v=V w (q) , eo). 
; Ig | 


Second Sound 


For an isotropic material, the velocity of second 
sound is (Niklasson (1970) ) 


-3 -3 
(vo + 2v, ) 


For hexagonal materials, the second sound veloci- 
ty is not so easily obtained. Meuller et al (1973) 


have assumed 
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which is an over-simplification. 
The phase velocity of second sound in hexagonal 
crystals is given by a theory in Appendix A2. 


The velocity Cy (8) for this theory is 
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where j is the mode (longitudinal, transverse), and 98 
is the orientation 6, (8) is the deviation of the pulse 
from travel perpendicular to the tranducer. 

These integrals have been arrived ne assuming a 
harmonic solid with linear dispersion using a Debye 


approximation. 
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Zz.3> © 6 Sound Velocity as a Function of Temperature from 
Thermodynamic Arguments 


Iet us. assume that it is possible to represent the 
temperature change in the sound velocity in terms of 
the adiabatic compressibility 


- 
a: i 


aes axa 


(2.51) 


where op is the density and K.(T) 


=i on V/OP) . the 


adiabatic compressibility. We can write 


2 2 pit’ eee , 
Veatr = ve (0) -= D a x07 (2.5-2) 
or 
KUT) (d9P/OV) 
ear = = (2.5-3) 
K_(T) (8P/9V) » ‘ 
Now 
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s Pia ny 
oP oV 
ac. = 4, # (3 fa (2.5-5) 
Pp V oT Vv oT Pp 


the derivation of the above relation for Cpy7 Cy is given 


in the appendix, hence 
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K. = K.{2 + a 2) a | 
Now 

(ir), ~~ (32), Gal, = lr], 
or 

K.(T) = (2) [1 + te - ¢| (2.5-6) 
hence 

K,.(T) > K, (T) ro, aLl. © 
equation (2.5-6) can also be written 

1 1 TV (aP)? 
Kay - K(ty + c [35 4 (2.5-6a) 


Now we have experimental data at several volumes for 
the right hand side, so let us write the above as 
1 1 cv. BBe* 


(ap ba. (2.5—7) 
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v“(T) -v os Faun K, (0) Gy (2) aT y 


If we calculate this value we find v2 (T)-v (0) is 


always negative and of the form 


v7(T) - v*(0) = ar? 


where n ~ 4 and A < 0 (Jarvis et al (1968)). 
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One can also proceed from (2.5-1) assuming a 


reduced equation of state 


S = S(T.¢(V)) = S(x) (2% 078) 


where $(V) is usually related to the Debye Theta by 
o(Vv) = — 


but we will not specify this. 


Now 


rl, > Gl, 
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Using (2.6-8) we have 
2 Oo a fe oe ox — oe o> (V) 
ds = 3 ax = 53/GR ars GH av] = 38/o var +11 VU ) 


32 


av] 
T 


$e 










6 pnrmunes ([-€.S) moat Beeoo71g ols aso Sad 


atetd2 Yo vod seups beoubst 
; eet 
(nea = ((Vje.T)e = @ 


' 


i 
f teiT syded of% of Begsior ylisuay el (V)@ sienw” Y 


4 


= (vie 
a® 


.2idd ytiosqe yon Iliw ow Sud 





‘ache sity va 





ss 








33 


oT aV 
hence 
35 36 
mel =). nO (VY) 
E Vv ax 
fi), = 28H), « he, 
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(38) ise 2 gn a 
aV , V oV . 
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Cp = c,, (T+ 9 (V) ) ' 
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Now if we combine (2.5-9) and (2.5-10) and inte- 
grate over T at constant V and if we recall $¢(V) and its 


derivatives are approximately temperature independent, then 
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This is identical to a result obtained by Edwards and 


Pandorf (1965) if ay/aV = 0 


2 
eo ee kee eck ee : 
Bey ok (0) AU; + 4-- Sy] = adu (2.5<12) 
If we make the association of ¢ = 1/6 , then y is just 


the Gruneisen y 


wm 1 OL 4 0.0835V (Ahlers, 1970) , 
hence 
2 2 a 
v (T) - v (0) = 5 [aulm> 0 for V < 24.2-cc/mole 
4 
~ AT for Pass BY 5, 


This result contradicts (2.5-7) which says A < 0. 

Our experimental results agree quite well with (2.5=7) 
(see section 4.3-2, figure 6). Thus a reduced equa- 
tion of state should only be used with extreme care for 


solid helium. 


2.6 Temperature Dependence of the Sound Velocity from 
Elasticity and Thermodynamic Arguments 


By definition, the adiabatic elastic constants 


are written 
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(2-6-2) 


(2.6-3) 


(2.6-4) 


(256-5) 


Here o,€ as subscripts mean holding all g,e 
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the thermal expansion tensor. 


for an isotropic crystal 
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is the dilatation. 
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(2.6-6) 


(2.6-7) 
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Now 


where xy is the thermal conductivity, or 


xver = Tr(S2) + cy (24) (yeas 


Now in general the equation of motion is 
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s component of the displacement 


where S; (x) denotes the i 
of the particles at point x. By substituting the above 


into (2.7-12) we get: 
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there are no losses to heat - thus truly transverse 


waves will propagate without heat losses. 


Lets try a plane wave solution to (2.7-10) of the form 
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but first as 
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make it numerically valid for hcp materials is difficult 
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2.7? Tsaund Velocity as a Function of Temperature When 
Coupled to the Localized Temperature Field 


The following is a brief description of a series 
of papers by Niklasson and Sjélander (1968), and 
Niklasson (1970) giving a theory of transport quanti- 
ties in anharmonic crystal using quantum statistical 
mechanics. Niklasson derives, and then solves a gen- 
eralized transport equation (2.6-1) similar to the 


Boltzmann-Peierls equation 
As (q,w) 
[2wa-M, (1) +M, (1') ]——>— n(w) (14n(w) ly; (q1w7 2,2) 


+ id, (a9) a5 = 2A, (4pw)n(w) [14m (w) 1H? (qui 92) <u (Q,2) > 


b2-ayi—1) 
where w is the energy of a phonon of momentum q; 

Q is the energy of a disturbance of wave vector Q; 

j is the mode (longitudinal or transverse) ; 

(1) and (1') stand for (q+%Q;w+s2) and (q-%Q? w-kQ) 
respectively; 

Aj (qiw) is the generalized spectral function; 

n(w) is the equilibrium occupation number; 

as (qu 7 Qn) is the deviation from equilibrium occupa- 
tion number; 

Hs (qui Qa) gives the coupling between thermal mo- 


tions and lattice deformation; 
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L, (2,2) is the generalized collision operator; 
<u(Q,2)> is the lattice deformation; 
M(t) is the diagonal part of the time-ordered 


equilibrium self energy. 


We will only be solving this equation in the 
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limit of low frequencies in the harmonic approximation. 


In the harmonic approximation, ws (q) are the harmonic 


phonon energies for mode j and momentum q. 


Aj(qrw) = F (6(wto;(q)) + 6 (w-w,; (q))} (2.7-2) 


nie) = (exp (hu/ET}— 1) * 


and if we specify 2 << ¢7H << 2 where t is a typical 


phonon relaxation time and w_ a typical phonon energy, 


a 
then 
2 
: = : 2.7-3 
st ee Ws (q) ( ) 
L5 (202) = £, (0,0) =. , the collision operator 
w% (q) 
: apne? ' ee 4 aa a, 

H (qwiQ2) = HOO (quiQa) = gr AL a ) Yap (a) 9, 


where v4 is related to the generalized microscopic 
Gruneisen constant ¥4 4) by 

as ji “ag 
Y; (a) = L. vg a anv 


(2.7-4) 
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where Uy is the dilatation and V is the volume. For 


a pure expansion or compression 


yi(q) = 4 Tr yII(q) . 


Now define the inner product of f and g by 


<f|g> = : | “4 | A; (q,w)n(w) (14n(w))£%g ; (2.7-5) 


Finally we will separate a function into q, w 
and j dependences by the projection operator onto q,wW 


as 


_ {|w><w a><qa = 
P= 2h ae i” l eee ae 6) 


then we can write 
a, = [a> = Jw>a ot |q>.a + (1-P) lar (2.4-73 


and rewriting (2.6-1) 


hi 


la> = (G6 + if) > |H>-u (2.7-la) 


where G is a flow term, [I a collision term. Now in our 
approximation 
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G=G =2 - +— cJ(q).9 
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where ¢J(q) = Ve w (4) the phase velocity . (2.7-8) 


If we now define 


pe (1-P) (G+ if) (1-P) 
and 


Rt = [(1-P) (G+ ir) (1-P)]7> 


(note: R+ does not exist when operating on |q> or |w> 
states; in fact Rt w> = BR 1q> = 0 = R|w> = R{q>), 


then we can formally solve 2.71 by 
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where 


Pej ie 0 
Oo 
relg> = 0 


u Abed 2 : 
thus i represents resistive (momentum non-conservation) 
scattering. If we further realize [ >> ee in our region 


(collision dominated regime) 


ea ek a if 
Oo 


reed. = - 

R [(1-P)T. (1-P)] as “To Pr Pra iG PY e<<T, 
(2.7-106) 

hence for cubic or isotropic crystals (hcp is dealt with 


in the appendix) 
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are different averages over the 


nitude as -¢ the total relaxation time. Cir is the 


velocity of second sound as usually measured and a 
is the umklapp relaxation times. These constants are 
calculated for the isotropic and hcp cases for He? in 


the appendix. 
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then we can write to first order in ane the 
generalized heat conductivity 
OO 
t+ — - i 1% 
Gey (0,0) .0 = eters os i elie ae (Dea 3 
i et eee alt 
tat Are diozi.OF7 thy(Q) 
i i 
u 
which in the static case becomes 
TOLLE LAs Gaqugtiog (2.7-14) 
O Veciieewus 5 ; ; 


Notice there is no Poiseuille Flow as we have ignored 
multiple scattering effects, The second sound velo- 


; : OO ' 
SLcty per To << ths 


A =3 
<c] (q)w. (a) |a,> ” a 
c¢.7.2 =—5 = for isotropic case . 


eis = 120:) 


Finally, the result we need, is the propagation of 
first sound with coupling to the localized temperature 


for the isotropic case. 
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This theory essentially presumes that around 
some localized distortion of the lattice, will appear 
a local distortion of the phonon spectrum. These 
phonons will be travelling in general with some 
group drift velocity v and decaying by various modes 
to equilibrium. 

The drift velocity in general is characterized 


by the following equation 


; 1 2 
bat Hel v (04m) + Co} Ut". 60 0.v, (0,2) 
a Oe II By6 aBbyd-y dé B 


35 ; 
am <W | W> Pa L iss 
= (eggs) [> CyQ,7 (2,2) 


+ ic, y2 dT ay sey 2g <Ug (21%) >] 
By6 

where the term in 1/t,, is resistive losses, the term 
in t" is diffusive losses, the term in T is coherent 
heat flow losses, and the term in t‘ is losses due to 
coupling between the drift velocity and the lattice 
deformation. 

The generalized heat conductivity is 
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If we use equations (2.7-12) 
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sound; 
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Hence propagation is enhanced due to coupling to second 


sound. 
ometor . o< Qty ven 


Hae (et 40) em 22 (0) } ; 


(2.7-19) 
J {r-O+ 7" (Q)} 


meme. 20C)) ee ef). 1 > then N. (Q) is usually negative. 

For se = 0; t"(Q) = 2tT'(Q) produces isothermal pro- 
pagation (usually 2s < 0 and we never have isothermal 
propagation). (Note: for pure transverse modes Q. =), 


hence Ns = 0 and propagation is always isothermal or 


s +eT 
always adiabatic as ci 7 Cj)- 


2.8 ‘Interference at the Transducer Due to Misalignment 


When the mirror and tranducer are improperly 
aligned, the reflected pulse will strike the transducer 
at an angle. The arrival time of the pulse is, therefore, 
no longer sharp. Also, interference across the transducer 
will affect the amplitude of the received pulse. 


Assume an input pulse f(t) where 


E(t) S"Arsin Le CG 25h ue 


0 otherwise. 
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Its fourier transform F(w) 
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travels with velocity V, the arrival time of the we pulse 


= , where w is the frequency of the pulse. If 
our tranducer is a square of size z, we must modify our 


inverse transform to be 
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This time change for the same parameters introduces 





an error in the direction of shortes times of 56° or 


.04% for two consecutive pulses. 
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CHAPTER 3 


EXPERIMENT 


eee LOM Temperature Apparatus 


Seb—] Cryostat 


Helium is the only element which does not solidify 
under its own vapour pressure. The phase diagram (fig. 8) 
shows that a minimum pressure of 25 bar is needed to 
solidafy Het. As we wished to study hcp He’ up to 
pressures of 150 bar, the cryostat was constructed to 
handle He4 gas, liquid, and solid at pressures up to 
250 bar (a safety margin of 100 bar). 

Whereas growth at constant volume tends to 
produce polycrystalline samples, the crystals for this 
experiment were grown using a constant pressure technique. 
The bottom of the pressure cell was connected by a copper 
strip to the He? fluid container. The pressure cell was 
thus cooled from the bottom. The temperature at the top 
of the pressure cell was kept slightly above the melting 
temperature (20 to 50 mk) during crystal growth by a 
322/£t manganin wire automatically regulated heater 
running inside the pressure capillary. This heater 
ensured that no solid could form in the capillary and 
block it, ensuring constant pressure growth. A thermom- 
eter on the top outside of the pressure cell allowed 


manual and automatic monitoring of the temperature. 
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Figure 9 illustrates the standard He* cryostat 
arrangement while figure 10 shows the pumping system. 

The pressure regulator on the He4 pot line (Walker (1959) ) 
maintains the temperature variations in che pot toa 

few mk between 1.4k and 4k. By placing a calibrated 
needle value on the ballast volume (Ackerman (1967)), 

one can cool the bottom of the pressure cell as eri 

as one wishes, facilitating automatic regulation of 
crystat growth. 

One filling of the He* pot lasts for about 16 
hours at 2k. However, the growth period can be extended 
by refilling the-pot through the fill needle valve while 
pumping on the fe pot. Although extreme care must. be 
used, the growth period may be extended indefinitely. 

A rough estimate of the temperature in the pot can be 
obtained by using a vapor pressure reading on a bourdon 
gauge (model FA 145, Wallace and Tierman, Belleville, 
N.J.) in the pump line. 

Onerritiing of the He? pot lasted from 30-60 
min. depending on the heat input. The lowest temperature 


achievable when the ultrasonics were on was about 0.7k. 


3.1-2 The High Pressure Cell 
Originally, these experiments were aimed at 
determining the temperature dependence of the ultrasonic 


attenuation. This consideration affected the design of 
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FIGURE 9 


CRYOSTAT 


Copper Can (22" x 3" diameter) 

Needle Value 

Thermometer-Heater (1002 Allen-Bradley) 

High Pressure Capillary (3/32" OD x 1/32" ID) 
He? Pot (4 1/2" x 2 3/4" diameter) 

He? Pot (1" x 1" diameter) 


High Pressure Cell (see figure 11) 


Thermometer (cryocal 2443, 
type CR 250 - 1.5 - 40 - He?) 


Vapor Pressure Capillary (1/32" OD) 
Capillary is Anchored to He# Pot 
He? Liquid 

He? Liquid 

Thermometer (1002 Allen-Bradley) 
Copper Bar (6" x 1 1/4" x 1/4") 


Vapor Pressure Bulb (2" x 1/2" x 3/8") 
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FIGURE 10 
PUMPING AND PRESSURE STABILIZATION SYSTEM 


A He? Storage 

2. Walker Pressure Regulator 
3. Ballast Volume 

4. Needle Value 


5. Cold Cathode Guage (H.S. Martin & Sons, 
Evanston, I11l.) 


6. Diffusion Pump (Edwards High Vacuum, 350 Watts) 
7. Exchange Gas Storage (He4) 

8. He? Rotary Pump (sealed) 

a. He4 Pot Rotary Pump 
10. He’ Exchange and Backing Pump 

R to recovery system 


V Vent to atmosphere 
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the pressure cell. A single transducer and mirror 
symmetry were chosen in that it allows observation of 
many echoes. The mirror had a diameter three times the 
mirror transducer distance of about 1 cm allowing 
deviations of the beam path from perpendicular of up to 
A595, Typical velastic.constants. (for example those of 
Crepeau et al (1970)) predict deviations between ray and 
wave normal of up to 20° for longitudinal waves and 35° 
for transverse waves. Thus velocities in crystals of 
all orientations were observable. The distance from 
mirror to transducer was chosen to be as small as 
possible while minimizing the effect of ringing by the 
transducer. 

Figure 11 shows the high pressure cell and the 
tranducer mounting. At first the transducer assembly 
was hung freely by three 2-inch long brass bolts from 
a thread in plate. Alignment was set at room temperature 
by examining the echo pattern in xylene and acetone. 
This system did not maintain the alignment to 4k. 

The system of suspension was therefore changed 
so that the suspending plate was no longer threaded 
in but rather held between a ledge and a split ring 
circlip tightly by six set screws. . The transducer 
assembly was suspended from this plate by a set of six 
push-pull bolts. Three bolts pulled the assembly toward 
the suspension plate while three bolts pressed it away 


from the suspension plate. This system maintained align- 
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High Pressure Cell 


FIGURE 11 
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ment well to 4k, however, repeated cyclings to 4k and 
room temperature gradually brought the transducer and 
mirror out of alignment. 

it insiaimtexestimes to'vpormmbtbwout thatkextneme 
tightening of the push-pull bolts caused distortion in 
the holding plate which gave rise to interference patterns 
in the ultrasonic echo pattern. 

The mirror was a copper plate 3cm in diameter 
(although only 2.8 cm was usable by the transducer) and 
one-quarter inch thick, polished flat to within 3 fringes 
of Na light. A range of polishing powders was used 
including various. grades of grit papers, alumdum powder, 
emery powder and rouge powder. After the mirror was 
threaded into the bottom of the pressure cell and soft 
soldered it was only flat to 10 fringes of Na light 
across its usable surface. This corresponds to roughly 

% of a wavelength of 5MHz sound in a crystal at 120 bar. 

The bottom of the cell (the copper mirror) was 
attached by thread and nut to a 12" x 2" x 4%" copper bar 
bent into a "U" and fastened to the base of the He? DOU. 

To the top of the pressure cell was attached the 
pressure capillary. The 322/ft wire was grounded inside 
the top of the pressure cell. 

The success rate for quality crystals in this 
cell was one in four at 120 bar and eva in three at 


85 bar. 
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3.1-3 Crystal Growth 


My experience with single crystal growth at 120 
bar and 85 bar appears to be slightly different than 
growth at lower pressures (for example Crepeau et al 
(1970), Wanner and Franck (1970) and Vignos and Fairbank 
(1966)). Growth from the superfluid phase is only possible 
at pressures in the neighbourhood of 25 bar. 

The growth procedure I have followed is outlined 
below. The pressure cell was brought to temperature 
equilibrium by letting it stabilize for about two hours 
with the sound system operating at a temperature about 
20 mk above the melting point. The automatic heater 
was then set to maintain this temperature while the low 
flow calibrated needle valve (Ackerman (1967)) in 
conjunction with the pressure regulator (Walker (1959)) 
lowered the temperature of the He4 pot. This pot was 
connected by He? fluid and the copper "U" bar to the 
bottom of the pressure cell. Thus the bottom of the 
cell was slowly cooled while the top of the cell was 
kept above the melting temperature for a period of 
12-24 hours. The longer the growth time, in general, 
the better quality the crystal. This final ‘temperature 
gradient was then held for another 24 hours to allow 
crystal defects not in thermal equilibrium (dislocations, 
for example) to move to the surface ewe the portion 
of the crystal used for ultrasonic measurements. The 


temperature at the top of the cell was then lowered by 
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slowly changing the setting on the automatic heater untid 
the whole cell was filled with solid. If the crystal 

was still low in echoes further annealing would be done 
(a further 24 hours) followed by cycling down to 1.3k 

and then back close to the melting point. If the crystal 
Wasestilltof insufficient »«quality, it was melted. 
Crystals of over 50 echoes were kept. Only one crystal 
with less than 50 echoes was studied, "D3', as it had 

20 echoes of very good pulse shape. This crystal how- 
ever, showed a knee at very high temperatures (2.6k) which 
may be due to the fact it was not good quality (see 
sections ?4eandm5)< 

The rejection rate at 120 bar was three crystals 
out of four, while at 85 bar was one out of three. From 
my experience it was easier to grow crystals at lower 
pressures. 

Although no direct evidence was available to 
determine whether the crystals were single or not, there 
was indirect evidence available. At 85 bar we measured 
sound velocities that varied up to 15%. This high 
anisotropy is identical to that predicted by Gilles et al 
(1968) and found experimentally at lower pressures by 
Creapeau et al (1970), Wanner and Franck (1970) and 
Greywall (1970). Also it is unlikely that crystals with 
low angle grain boundaries or small crystallites could 
produce over 50 ultrasonic echoes (attenuations of the 


order of .l om Ehs For these reasons we believe that 
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our crystals were single. 

The ultrasonics were left running during growth 
to facilitate monitoring. One observes a decay in 
amplitude of liquid echoes until they almost completely 
disappear. At this stage, probably, large numbers of 
crystallites have formed and covered the mirror giving 
random reflection. Then, slowly, solid echoes begin 
to grow (displaced to shorter times). Annealing 


improves both echoe quantity and quality. 


3.2 Temperature Measurement and Control 

The temperature at the bottom of the pressure 
cell was determined by an Allan Bradley resistor (100 %) 
whose resistance was measured by an Oxford instrument 
resistance thermometer bridge which was rigged for two 
and three terminal resistance measurements. This allows 
one to compute the four terminals resistance assuming 
the resistance to all four leads is identical. The 
resistor was prepared by grinding the outer coating down 
to the carbon centre and then covering with a thin coat 
of GE varnish (#7031) and baked at 90°C for two hours. 
The resistor was then inserted into a tightly fitted 
machined hole in a copper block. The lead wires were 
wound several times around the end of the block, then 
wound around the copper "U" bar, then wound several 
times around the He? pot. Finally the leads were anchored 


at 4k and then taken through a vacuum feed - through to 






. ay 


a 
my : 


> 


-sipnis sxew eleteyto 


o 
§ 








WOIp pritol patanux tiel susw epinoes1siu oA? 


. rm : 
vVrsedo ent -PAL TOA inom st3e7 ii tos? oF 
aes 


‘4 ——— e * aor s —_ : 
8 yous lijnw esorse biupii to shoutilen 


; 2si .yldsd 21g ,epase wine 2A .aesqqsekh 
& 







sieved boas beorot evad. estiffea 


-(29mMiz? zwedsgiona oF beoxs sige fb) worp. oO: 


- : 4 —— , ‘ 
DMS yIisnsup sodoe fijod 







zsvoxqmt 












7, 
. . 
_ 
~ BbaB InometTwascaM ors syeqneT Sat 
fod sdj 3s IIIS teqms? ad?T 
~ , a 
. 
. oe ff t yr ie > - zg 
; b|7a nal! 1% YG benterreteb esw 99 


; 2x0 m6 yd bowesem esw eonsteizes1 seonw 
10% bap Saw doidw 651d wseademonyroeriz en 
£ : a 


isivesem sonstdeiear ILaniozes ontdt 

16 sonstelest elsoimas? Wo? aris. stuqmoo oF > 

7, 

sAT s0itinebi ei absel.2zwot Ils of oo1sseisex’ 6 


mwob paisso> 1e3u0 efi pribnize yd bexaqerq asw 2092 





it alas ptats aa ealteves seit brs e13ne8D nod3s5 


ma 





ie 


* ows °C . >*0e. j r92 - 
panes: | - ce ee - 's 


70 


room temperature at the top of the cryostat and connected 
to the bridge. 
The new resistor at the top of the cell was also 
a 100 2 Allen-Bradley prepared in the same way and 
attached to the top outside of the pressure cell by 
bolt and grease. It was not anchored at the copper "U" 
bar as it was used during crystal growth, but was anchored 
at the He’ pot and at 4k. 
The thermometry was later changed at the bottom 
of the pressure cell using a cryocal, Inc. calibrated 
resistor CR 250-1.5-4d-He~ S/N 2443 Germanium. This 
resistor was coated with a mixture of vacuum grease and 
copper filings and inserted into the "U" bar at the 
bottom of the pressure cell. The resistors were calibrated 
against the 1962 He* scale and the 1958 He4 scale of 
temperatures by a vapour pressure cell bolted to the 
bottom of the —"U" bar. 
AW the top of the He4 pot two resistors were 
used. One was a roughly calibrated two terminal 
thermometer (calibrated using the boiling point of He* 
and the }i point) and one was a simple 100 2 heater. 
A 302/ft manganin wire was placed inside the 
high pressure capillary and grounded to the top inside 
of the pressure cell. This heater was used in conjunction 
with the thermometer at the top outside of the pressure 
cell in the following manner. A resistance R correspond- 


ing to the desired temperature is set in the temperature 
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Temperature regulator schematic 


FIGURE 12 
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Fequbator~;(21dure 12) “liethe cell, is “too: cold, anve.m.f. 
is set up across the amplifier (105) which turns on the 
transistor allowing current to flow through the manganin 
wire. When the top of the cell reaches the desired tem- 
perature the e.m.f. goes to zero and the transistor shuts 
off. By choosing the time constant properly, the cell 
can be continuously cooled while maintaining the temper- 
ature at the top of the pressure cell constant. The 
"Manual-auto" mode switch allows the proper settings to 
be worked out, without damaging the crystals. 

I would like to thank Len Vienneau for helpful 
discussions which led to the use of large temperature 
gradients during annealing. 

The calibration technique for the lower thermo- 
meter is accurate to 10 mk throughout the range of 


measurements. 


She: High Pressure Gas Handling Equipment 


The purity of the helium gas used is important 
for many types of measurement on He4 systems. Purifying 
by super-leak is undoubtedly the best method, although 
second sound has been observed in He? crystals prepared 
from gas evaporated from liquid helium. For measuring 
sound velocities, such high purity is probably unnecessary. 

Figure 13 illustrates our gas handling procedure. 
The system was first evacuated, then commercial grade 


He” is fed into an activated charcoal trap cooled by 
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liquid nitrogen to remove water vapor and other materials 
that will freeze at liquid nitrogen temperatures. The 
helium was then passed through a pressure trasmitter 
(necessary to achieve pressures greater than 180 bar). 

The gas was then fed to a nitrogen trap and 
allowed to reach thermal equilibrium. Finally the gas 
was Slowly fed through a liquid He* trapecill tagain 
thermal equilibrium was reached. The gas was then stored 
with the residue being pumped away after warming. 

The ideal limit of purity of gas prepared in this 
way is the natural abundance of He? in He’ (about one 
part in 10°). The gas prepared in this way was found 
satisfactory for the growth of single crystals and for 


the measurement of sound velocities. 


3.4 Ultrasonics 


3.4-1l Transducer Assembly 
Tne etalon ic a longitudinally cut )(xX-cut) 


quartz transducer flat to one part in 106 


coated with 
a thin film (100-1000 A°) of evaporated aluminum. The 
fundamental frequency for the etalon was chosen to be 
5 MHz. 

The etalon was backed by a brass piston polished 
flat by hand to 4 a fringe of sodium light over its 


entire backing surface (1.2 cm). The backing piston is 


about 2.5 cm long and is inserted in a teflon sleeve. 
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The diameter of the piston and sleeve is the same as 
trac or the etalon (1.2 cm). «The piston’ is, spring loaded 
into a large brass cylinder (2.5 cm in diameter). The 
brass cylinder and piston were hand polished in a special 
assembly over 5 cm in diameter, so that both pieces could 
be polished simultaneously without rounding at the edges. 
The two assemblies were made flat to about 2 fringes of 
Na light by polishing with various alumdun, emery and 
rouge powders. 

The assembly was fronted with a disc flange 
2.5 cm in diameter and 1 cm interior hole,» 2°) mm thick. 
This brass flange was placed against the transducer to 
hold it against the spring. It was bolted by three bolts 
to the outer cylinder. This flange was also polished, 
Mowever, it was difficult to get it flat near the inner 
edge. | Thus, Lt was flat to one fringe of Na light 
except near the inner circle where five closely spaced 
fringes could be seen. Next the inside, polished part 
of the flange was slotted six times with a file to 
produce radial slots about 0.01 cm in depth. The 
purpose of the slots was to prevent liquid helium from 
being trapped close to the transducer. If this happens 
the liquid may freeze later than the rest of the crystal 
introducing possible defects close to the transducer, 
or even not freezing at all. 

This assembly was then mounted by push-pull screws 


to another assembly described in section 3.1-2 and is 
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idlustrated in figure 11. 


8.4.-2(% Electronics 

A block diagram of the electronics is illustrated 
in figure 14. The pulse originates in an Arenberg 
pulse generator (Model PG-650C mod IJD ADD5D SN704) which 
is pulsed by a Hewlett Packard (222A) pulse generator. 
The radio frequency signal was then matched by a stub 
matching network through a tee. One branch of the tee 
travels through the high pressure capillary and connects 
to the back of the piston in the etalon assembly (section 
3.4-1). The other branch of the tee takes the pulse 
through an Arenberg tunable preamplifies (PA 620 SN 312B) 
and then through an Arenberg Wide Band Amplifier (WA 600-D 
SN 298). The final result is displayed on a Tektronix 
Osetbloscope (R 556) or on a chart recorder (Hewlett 
Packard Mosely 580M). 

The initial pubse-was, from 20 to 200 Vpp, at 
5MHz..with.a.pulse widthefrom 0.58to~5 ys and a» repetition 
rate of 100 per second. The pulse generator and the main 
trigger of the scope were triggered by the the same 
pulse from a Tektronisc Time Mark Generator (type 184 
SN 5655). The time mark generator also produces a 
second trigger pulse coherent to the first pulse to 
one part in yes It was found that the internal delay 
trigger and the internal normal trigger had variable 
AtiieeeO Lk tio us an hour. External triggering eliminates 


this source of error. 
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Figure 14 is a block diagram illustrating the 
electronic configuration. The coherence method of 
measuring the velocity employs the fact that a change 
in velocity changes the start time of each pulse. As 
the pulses are all starting with the same phase, the 
measurement of the change in phase at a given time from 
the first pulse will be an indication of the change in 
velocity. The no pulse will record a change in time 
of ang where d is the distance from transducer to 
mirror and V is the velocity. .The larger the n, 
the larger the velocity change. Sensitivity is therefore 
enhanced by using later echoes, unfortunately attenuation 
limits the echo so chosen to about n= 15 to 25. The 
time chosen for watching phase shifts must be very stable 
relative to the initial pulse time. This accounts for 
the care taken to provide this stability. 

One wishes to transfer the maximum electronic 
energy at the generator to sound motion in the crystal. 
The acoustic mismatch between solid He* and quartz 
is quite high (about a factor 100). The transducer 
is equivalent to a 1 M2 resistor ina serreswieh 
a 20 pf capacitor. The high frequency line in the 
high pressure tube has a transmission line impedance of 
100 2. The output impedance of the pulse oscillator 
is 50 2. The input impedance of the amplifiers is 
90 2. To match the cable in the high pressure tube to 


the transducer directly would be very difficult as well | 
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as only being exact at one frequency. The method which 
gives good matching without being terribly sensitive to 
frequency is the stub matching network being composed 

of three lengths of cable like the letter I with 

input and output at the two vertices. The theory needs 
the length of cable between source and receiver to be 
comparable to the wavelength of the cable. For our case 
the wavelength was about 40:m and the cable length 
about 10 m. Fach of the two hanging cables were terminated 
with a variable inductor in parallel with a variable 
capacitor. This allowed the network to be matched with- 
out cutting the cable. The solution for the input 
impedance . Zs of the network for a load Zs and cable 
length ‘a' and wave vector k is 


Jig ee 
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Te e 2ika) 
where Zo is the impedance of the cable. Single stub 
matching is a 'T' type network, while double stub 
matching is a 'Il' type network. Although the values 
can be worked out in theory, the various impedances are 
not well enough known to be calculated in practice. Thus 
we used a double stub matching network terminated in 
variable capacitors and inductors. 

It was found that the matching giving the largest 
echoes would change with temperature. Although the 


velocity measurements were not affected by this, the 
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attenuation measurements were made more difficult. 


3.4-3 Method 

The first step in measuring velocities was to 
calibrate the distance travelled. This was done by 
finding the time of flight in our pressure cell at 
various pressures in equilibrium with exchange gas to 
the Hea bath (4.15 k). This result was then compared 
at several pressures to known velocities in the liquid. 
The results are discussed in 4.1. 

The pulse echo method was used to calibrate the 
distance as well as to establish the absolute value of 
the velocity in the solid. The transducer was excited 
with a pulse from the pulse generator (20 to 200 Vpp 
for 5 to 5 ws with a repetition rate of 100 pulses per 
second). The transducer then mechanically vibrates at 
5MHz and sends a sound pulse through the crystal which 
is reflected at the mirror and returned to the transducer. 
Part of this pulse was then reflected and part transmitted 
to the recording equipment. If one then rectifies and 
integrates the signal using a gate of about 1 us, 
one gets a series of pulses out as in figure 16. The 
time delay between each pulse was measured. To calculate 
the average time of flight, a least squares fit to 10 
or 20 echoes yielded most accurate results (this method 
tends to nullify some systematic errors compared to a 


simple mean, for example). Ringing of the transducer can 


>a 
=< 
















») 
| 7 
; ‘Oon 4 inamreyTrens : ne it 
renem noijeunsiie | 
| 
- r 
Ol i3oM t- BLE 
j 59m mi gad ti2 oT 
bea flavs7s scunes > SAZ 
e| >t ' - r+ 


. in UB OL svuolisy 


; ‘7 ~ta €f.6 fijsd oH eZ 


2 siq iszevoe 
| : : 
ilivuesy sf 
A... } f 
< AiG ernT 


4 

“OW 86 Sonsiet 
‘ Zs + } P a 

i. - SAT AL vo (20lsv, odd 
: 7 i? mor? selug s dtiw 
ad 

? 


ey ¢€ og @ 107 








t vbens33 ocT . (Bnoss 
: a 
¢ 2@ & #Dnee Bana st Ma 
% 
; Pe 7 itim sf ¥en betosltex aj 
si g 


q (e135 20. 3} 


bas seltisne:x nodd ano I inomyl ps 1 
* a mT PS p card ta 3 
*] 


sirpie Mia ial i 


> i 










81 


cause distortion in the first few echoes. For this 
reason, it is better to take measurements on later echoes. 
We chose echoes around the twentieth for almost all 
measurements. 

Once the absolute velocity has been determined, 
the phase comparison method can be used to determine the 
change in velocity. For this method a pulse with good 
shape was chosen and one with which the radio frequency 
signal was undistorted. The amplifier was adjusted to 
transmit a non-recfified signal and then part of the 
echo maximum was chosen to be centred on the scope. 
Because of the external triggering on the scope (every 
10 us), only one part of each peak is accessible. Our 
appropriate scale was now chosen to ensure linearity of 
the scope over the measurement region, as well as 
convenience for measurement. It was found that an 'A' 
scale setting of 100us/cm and a 'B' scale reading of 
(0.05 - O.lus/em) gave best results. Normally, the delay 
was in the 500 us region (the twentieth echo). At our 


7 seconds or two to four 


frequency a period was 2 x Loy 
cm. on our scale. It was found that the signal could 
be read to 1/50 of a wave length. As the temperature 
was changed, the pulse would move to lower time of 

flights and the change in phase was recorded. As the 
change in phase is proportional to echo number, this 


had to be taken into account in the calculations. The 


appendix contains computer programs used in the numerical 
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analysis. 


3.4-4 Error Evaluation 

Errors in the pulse-echo technique came mainly 
from two sources: variation in pulse shape and the 
linearity in the delay trigger. The "A' scale was 
maintained on one scale to prevent errors by scale 
changes, A typical-time of flight was 25 us, so to 
see 40 pulses one had to be on the 100us/cm scale. The 
linearity in the delay trigger is quoted as_0.5%. Errors 
in the pulse shape were much more difficult to analyse. 
That part of the pulse shape error that is systematic 
can be removed by taking a least squares fit to several 
pulses. However, random contributions, as well as 
systematic errors of frequency of the order of the 
sampling time are more difficult to evaluate. Hope- 
fully ringing effects were small as we were taking our 
pulses at least 15 or so from the original pulse. There 


is a possibility of the pian 


transverse branch being 
propagated along with the original longitudinal pulse. 
This is very difficult to remove. Random scatter 

from the sides of the bomb as well as the effects of 
misalignment and diffraction also will contribute to 
ere The best indication of error, is therefore, the 
standard deviation of the least square fit to the echoes - 
normally about 3%. As the error in the distance calibra- 


tion is approximately 1%, the total errors were probably 
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The limitations of the phase detection 
technique were the flutter time of the oscilloscope 
triggers, the period of the radio frequency oscillation 
and the distortion to the signal. 


For a frequency of 5 MHz, the period is only 


Ze x fare seconds. If one can read the change in phase 


to 1/50 of a wave length,one can achieve a precision of 


4 x ae seconds. Compared to our mean flight time of 


ee x 1077 seconds we have an accuracy of 1.6 x hos: 


The flutter of the delay trigger is 2 x LOT 


an error Of 1.7: x Omer At low amplitudes (at low 


which gives 


temperatures the attenuation was high) the signal would 
be slightly distorted, perhaps doubling the quoted error. 
Some systematic errors are not important here, as we are 
measuring only the relative change (for example, the 
error in distance calibration does not contribute to 
this error). Greater accuracy could be obtained by 
using higher frequencies, but this would still be 


5 


limited by scope flutter (5 x 10 ~). ‘The phase 


superposition technique is even better, allowing 


5 6 (mMcskimin (1961)). This 


accuracies of 10 ~ to 10. 
method can be extended by integrating over the spectrum 
and using two pulses, one slightly delayed (R. Wanner - 
private communication). Using more pulses is possible 
also. J. Holder (1970) has devised a method using a 


large number of pulses properly delayed and integrated 


cc) 
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to accuracy of Mohseee All these methods are presently 


being investigated, although the latter would cause 


overheating in the sample. 
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CHAPTER 4 


RESULTS 


4.1 Sound Velocity in Fluid He? as a Function of Pressure 


In order to calibrate the distance between the 
transducer and the mirror, the sound velocity was measured 
from 25 bar to 135 bar in the fluid at 4.15k. This 
velocity was then compared to the corrected values of 
Vignos and Fairbank (1966) from 1 bar to 50 bar in the 
region of overlap. 

The pulse echo method was used to measure these 
velocities. A sound pulse was sent out from the 
transducer to the mirror which was reflected back. The 
electronics detected the echo as well as any further 
echoes and the echo train was then displayed on an 
oscilloscope (Tektronix R 556). For a more detailed 
description of the technique see section 3.4-3. 

The high pressure cell was maintained at 
constant temperature by thermal contact via He* 
exchange gas to a liquid helium bath at 4.15k. The 
data of Vignos and Fairbank (1966) was extrapolated to 
4.15k using a quadratic fit to their data at from .5k 
to 4k. Then a correction of + 1% was added to 
account for a systematic error noted by Wanner (thesis, 
1970) and Abraham et al (1970). Using these values in 
the region of overlap, a value for the distance was 
determined. 
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Figure 15 shows the resulting calculated 
velocities plotted with the results of Vignos and 
Fairbank vs the square root of the pressure. 

Notice that the relationship is roughly a 
straight line. 

Table 1 lists the results of the present series 
of experiments. The best fit of the velocity V(m/s) 


to the pressure P (bar) is 


= * 
Vi w= Ap + A,P + AP 
Ay am 130 16635 
Ay = 447.2186 
A, = -0.472669 


A.2 The Echo Envelope 


The echo pattern in an absorptive medium with 
properly adjusted mirror and transducer should be 
decaying exponential. The decay rate is an indication 
of the attenuation in the material. 

Figure 16 (top) shows the echo pattern for 
acetone in our pressure cell at room temperature after 
alignment at one atmosphere pressure with a separation 
between transducer and mirror of 1.46 cm. The velocity 
at these temperatures is about 1200 m/s. and the pulse 


can be seen for about 3ms. 


87 



















JI6ivoleas pritives: sis ewore 21 empl 

~- s , ru - a<6 T wees 7 ; 1 
| 2tinest edt dttw betsolq esitioolsy 
’ “7 aa « 


‘USeStg osfit Jo Joou srsupe oft av Ansdti 
dtist ) 


a de 


vidpvor at qifenottsisa of tans onksor 


.snil Jnpis3izea 
. 
asx sriz ezJe ri { > fdeT 7 


2 \) oisv ef? 10 34% Jeed sat .einemixsqxe 36 


ei (tad) ¢ sxwee “od 
) siveesig siz J 


+ q a whe fi = Vv 

ef 3 + ¥ 
i 

+> Ps = 


r ~ 
of 4 . \ * (A 
ee . 
a> a 
\ 
; 
scoleyvna odes off -S.8 
— lt. een tee a 7 
: @ 
i Gioceas MK AL aArs3tsq odos sent 


~ 
Oo bivor is ube 525 bos » or, im bateuf bs yixego: q 
a oe 





no ‘1 te riba re oT Steyr Veneb' ott ‘us 7 a 7 ‘ a» a ~~ 
(290D's -fsitnenogxs palysoeb 

¢ 2 f . ~~. b 

«Isizweism offs ai notsisunstts\ it. to 





402, mr9335q once » 


7 & 7 


2 eworle (qos) 21 om 


* es ee ; a4 _ 


elt 














Pressure (Bar) 


L306 


EOO:.0 


TABLE 1 


Velocity. (m/s) 





G5 325 
50137 
Sys ieks Sis) 
Sd. G7. 
439.74 
439058 
415.39 
413.97 
SOF aa 1 
396297 


hs U2 


+ 


1% 


88 
























#6 
=a} 
i AJ@AT 
Si + (e\m) ytiselsyv (168) swesert en 
ec. ana a.0tL 
Te. 182 0.00L peer ton 
2£.022 ¥.ee : ae 
t3.£12 . rer 
BT. eee a.e 
82. een ©. 
ef.210 ee | 
re.efh a. _ 
1, Tee fae i 
Te. aee 0. ae = 
CT. ITE a.8s s 
d 1G 
cam 
: \is ae" | 
dl bot wi _ _ 


v 







OY oa Oe 





J er) ? / 
=e i Pa ous 


7) 
nad oe be bre oe " 






89 


119.6 BAR 





ea 


FIGURE 16 Echo Evelope: at upper for acetone at room 
temperature, and below for He solid at 3.7k. 
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Figure 16 (bottom) shows the echo pattern for 
hep He* at 3.7k and 120 bar for crystal D4. The maxima 
and minima correspond to interference as predicted in 
section 2.8. There is not complete extinction because 
of attenuation. Thus, even though the cell was properly 
aligned at room temperature, the cell was unaligned at 
eit Ke 

Error due to curvature of the mirror will only 
be of the order of seconds of arc whereas actual mis- 
alignment is of the order of minutes of arc. At present 
it is presumed that differential contraction causes the 


misalignment. 


4.3 Velocity as a Function of Temperature in Single 
Crystals of hcp He’ 

The velocity as a function of temperature in 
Single crystals of hcp He* has been measured at two 
volumes corresponding to melting pressures of about 
120 bar and 86 bar. 

The initial velocity was measured ae the pulse 
reflection method to an accuracy of about 3%. The change 
in velocity was then measured by the phase comparison 
technique to an accuracy of about 2 x ida These 
techniques, and the method of crystal growth are out- 
lined in sections 3.4-3 and 3.1-3. 


Figures 17 to 26 are graphs on 10 of the 12 


crystals measured. The relative velocity is plotted 
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against the absolute temperature. The relative velocity 
is the actual velocity divided by the velocity at the 
coldest temperature. 

The crystals are numbered such that a different 
letter implies a recyclying to room temperature. Thus 
crystal E2 would be the second crystal grown during 
a specific cool. down period. F2 would be the second 
crystal grown during the next cool down period. 

Most of the crystals grown show an abrupt change 
in slope in the temperature dependence of the velocity 
at about 2k for crystals at 120 bar and at 1k for 
crystals grown at 86 bar. We call this point the anomaly 
or knee. At temperatures higher than To (the temperature 
of the knee) the velocity roughly obeys a power law of 


the form 


Ves VCO) AT 


where n is about 4 (as predicted classically). 

At temperatures below T, the velocity is 
constant. Not all crystals grown show the anomaly (for 
example, crystal A3). However, the crystals with the 
most echoes appeared to have anomalies. 

Comparison of our results with various theories 
will be given in the next section. Complete data on our 
crystals are tabulated in the appendix. 

It was a possibility that the knee was caused 


by some defect in the measuring apparatus. That is, the 
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reason the velocity remained constant was because the 
temperature of the crystal was not changing. 

Firstly, it should be noted that the pulse height 
(that is, the attenuation) did not reach an extrema until 
-lk or .2k lower than the anomaly. This is evidence of 
the temperature changing after the anomaly. Secondly, 
the crystals grown at lower pressures have the anomaly 
shifted to about 1k. It thus appears that the cryostat 
does achieve temperatures in the crystal lower than the 
anomaly. 

It was also possible that our sound beam was 
heating the crystal to the point of making the temperature 
roughly constant. This heating, in general, would depend 
on the conductivity of the crystal and thus might be in 
a different spot for different melting pressures. 

The sound beam was 140 Vrms for 300 us per 
second and at a steady current of about 2 mA for a duty 
cycle of 80 wy watts. At 120 bar the thermal conductivity 
at 1k is about 10 watts/cm°k and at 2k is about 0.1 
| watts/cm°k. In our cell there would be created a temper- 
ature gradient of about 1 mk at 2k and 107? mk at 1k. 

To make sure heating from the ultrasonics was 
not a problem, an experiment was done with the duty 
cycle reduced by 100 (figure 27). The two curves are 


essentially identical. 
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Relative velocity as a function of tem- 
perature for crystal F3 at two different 
ultrasonic input powers. 
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This experiment was done on crystal F3 about 
5 days apart. Annealing might account for the slight 
difference between the two crystals at T > Tot 

A factor in favour of the effect being real is 
that A3 does not show the knee. Further there is large 


variation in the knee position as illustrated in table 2. 


4.4 Observations of the Attenuation 

To my knowledge, there is only one mention of 
the attenuation in an ultrasonic experiment (Vignos and 
Fairbank (1966)) in which an upper limit was placed on 
the attenuation \O€ jorder!}//3)to\:7 om7* at 10 MHz 
over their measurements in the solid from 1 to 150 bar. 

No accurate measurements of the attenuation 
could be made on our system due to several factors. 
The fact that the transducer and mirror were slightly 
misaligned caused disruptions in the echo envelope. 
It was also observed that two neighboring echoes would 
vary in size making the second larger than the first 
(or vice-versa) as a function of temperature. Also, a. 
general deterioration of pulse shape consistently 
occurred at low temperatures (but pulse shape was 
recoverable when the crystal was returned to higher 
temperatures). 

One of our crystals (crystal A8) had very 
‘small misalignment as well as a large number of echoes 


(140). If we blindly apply a least square fit of the 
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CRYSTAL A8 
120.8 BAR 
a7 kK 
ATTENUATION - 0.07 cm7! 





CRYSTAL A8 
120.8 BAR 
tek 
ATTENUATION - 0.56 cm”! 





FIGURE 28 Echo Pattern for crystal A8 at two 
temperatures. 
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form 
€(x) = Be. 


to the echo pattern we can get a rough idea of the 
attenuation and how it is changing. 

Crystal A8 was one of the first crystals in 
which a knee was observed. However, as this result was 
unexpected, proper care was not taken in the velocity 
measurement and systematic errors occurred. There 
was left a fine integrated pulse record, fortunately, 
and that was used in the present analysis. 

Table 3 shows the values of the attenuation at 
the two temperatures measured (figure 27). If we view 
these figures as an upper bound, the attenuation at 5 MHz 


1 


and 120 bar in solid He? is about 0.07 cm at 3.7k and 


Gase.cm 


at 1.3k. This apparent increase in attenuation 
was observed in all crystals except D5 (where a small 
peak may have been observed near 1k). 

In crystals showing a pronounced anomaly, the 
attenuation increased very rapidly in the knee region, 
and then increased more slowly to lowest temperatures. 


In all crystals the attenuation changed by about 


a factor of 3-7 from the melting point to .75k. 
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CHAPTER 5 


COMPARISON OF RESULTS TO THEORIES 


5.1 Sound Velocity in the Fluid 


In section 4.1 we showed experimental results 
giving a pressure "P" dependence for the velocity 


wy" of 


L 
WaéP* ; (5.1-1) 


In a corresponding states treatment we can 


write for He* 


wee Wa aa (14+ F*D*4+Gepe-4Hepe ry |} ; (51-2) 


where Pp* =P = o*P/e for Lennard-Jones Potential 


0 
of type 


yo eeede(ger) = (d7r) (513) 


For He4 it has been determined from second virial 


coefficient data (Hirschfelder et al (1954)) that 


ste k = 
Pp, = 87.4 =v = 89.1 BAR . 


_In the limit of P* = 0 we should get a value 


for le that is about 2 (normal for monotonic gases). 


2 
y = MWI(0) _ (147) 74 
Oo RT x “cl 


= 109 = 
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This particular treatment also assumes that P* 
He small. Incur case P* ~ 1 and as the value of i 
is not correct, this simple picture must be abandoned. 
De Boes and Lunbeck (1948) introduced corrections for 
quantum behavior by introducing an extra expansion in 


the quantum parameter |\* = st = 2.64 for He* 5 


ov me 


Because A* is greater than one, this expansion is 





invalid in our case. 


In a Similar fashion, all the simple theories 
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which I have examined dealing with dense fluids, includ- 


ing the above and the solid sphere available volume 
theories, give bad results in either the PVT relation 
or in the pressure dependence of the velocities. The 
theories in general predict higher powers in the 'P' 
dependence of the velocity (like p? in the Carnevale 
and Litovity (1955) theory) as well as much smaller 
changes in velocity with P than observed. 

From the above, I conclude that a proper under- 
standing of the velocity dependence in fluid He4 under 
pressure will depend on the development of a more exact 
theoretical treatment which takes into account quantum 


effects including short and long range correlations. 


5.2 Theoretical Treatment ef Velocity Results 


Belek. introdmetaon 


Figures 17 through 26 show the typical behavior 
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of the relative change in velocity. We will first try 
to explain those crystals showing the anomaly and then 
those which do not. The temperature where the abrupt 
change in the rate of change of the velocity (the knee) 
will be labelled To: This temperature, although not 
clearly defined due to discrepancies in the velocity in 
the region of the knee, is called To for convenience only. 
We will examine the data in two regions T.> Te. and 
es To: Above Toe the velocity appears to be adiabatic 
and reasonably classical. Where T § Ty we believe 
that coupling between the temperature field and the 
lattice field are taking place. One could alternately 
say that, in this region, temperature differences are 
being equalized by coherent propagation of heat waves, 
rather than diffusively, that is by second sound. 
Notice figure 29. The second sound region 


normally occurs fora T such that 


QT. ove BORE 8 QT, : 


The knee occurs roughly at QT 2 an anes cation 
that second sound propagation may be the key to under- 


standing the anomaly. 


De2.2 Velocity Change Wich, ii an 


Sections 2.5 and 2.6 give detailed theories for 
classical behavior of a solid. Briefly, I will mention 


their basic conclusions. 
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Classical theories all state that at high enough 
temperatures (or low enough frequencies) all sound motion 
will be adiabatic. Whenever a longitudinal wave is 
excited it produces compressions and rarefactions which 
in turn produce temperature gradients along the path of 
propagation. 

For example, a fluid with cross-section A to the 
normal of wave propagation, has, say, created a temper- 
ature gradient AT. Then the amount of heat Q, flowing 
across a cross-section in the time the wave travels })/2 
is 


Nei AT de SSR AT 
me ow 


Now the amount of heat required to equalize the temper- 


ature is Q.5 


Ar e: " AE COAST 
Smee 5 ay ae 
m m 


where»+k is the thermal conductivity; 
25 is the heat capacity at constant pressure; 
Ve is the molar volume; 
Vp -is the Debye velocity. 
veces 
ls kee Us < Th aa the propagation will be 
m 


adiabatic as opposed to isothermal. 


Now 
(Eien 23) 


or 
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Re. 
Cua pe Be ie 
er <a qT re) = uf e (522 4) 
Vp V 


From thermal conductivity data we know that 
Sy ich tem. tke tO LViSK ror SMH2 at? Bo* ear 
aid t.ok tO ls /k Lor SMHs ar ¥20' bar. 
Figure 28 shows mean free paths at 86 bar. Thus for 


T > 1.7k we should have adiabatic propagation. 


Cc 
Vann ee Va (T) +a(T)(- 1) wir) . (5225) 
V 


Writing the velocity in the form discussed in section 2 
we recall that V(T) is adiabatic for N(T) = 1 for all 
theories. To get an idea of the form of V(T) lets recall 


equation. (2.5-12) 





a a Vo di he Ns 
Few Mie es ave 


all velocity measurements are made at constant volume 


so that 
7 4 
U-U, = | CAT exrak for a Debye solid . 


If we assume that y is roughly constant in temperature 


we can write the adiabatic velocity as 


aie ks 4 


Further, for a Debye solid, the isothermal velocity 


._ > >> TH 
we 













S9ab yItvistoubaoo Isaryeiy mox 1 
8 +6 sHM@ yo! A#é@.1f of axAL.i ge { = a 
. y 


- m . 
=< OF i I AWS 
- 


6d 98 Js arijaq sex? asem eworle BS exupk® 


tq visedstbs sved Bbluode ow AT.I < | 
» ; =e 


.] La rf ; 4 ” => S 
i (L- $2) (re 4. (Ti Sv = Cntve Pie 
be 
a Ac ~ rt By: ea2u5i ‘5 miro} gr rr y2 tsoleay 4eri4 pnig 
I (T)4 Or Sf36d6i56 ef (TIV sede [issex sy 


: + 


- : > A - 
os oa i , Aa DS O 690Di ds J Sp ¢ T -eoit os 47 
» ee 


ia 
| (8E+8.8) aokssupe 


° 
i . mn? « 4 J 
= — ty a -_ = 
Ve “yt BUY) = yay ~ ee 
- . 5 { P an 
G 2 


. Biloe syded s x02 T » Th2 
: © 






(| ne 
 ~SIsIB 1e¢ 






115 
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20.72 cc/mole 
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FIGURE 30 Experimental for (c,/e -1) according to 


Jarvis et al (1968) 
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Now, aS was shown in section 2.5, helium does 
not appear to obey a reduced equation of state 
for the velocities. However, it was determined 


that in the region T > Toe 
n 
Vip (T) a A + Br 


Where 1 Uvaries from 3.5 to 4.5. 


If, as in the Debye model, we give Vor Var the 


c 
same temperature dependence as a we get 
V 
Zits n 
Vin = A + BT 
Cc 
Pe te Daa P it n 
Nab Ca Ly const: =)c.'+ DT 


as the temperature dependent terms are small. 
c 
Now figure 30 illustrates (= - 1) according 
V 


to Jarvis et al (1968). One can fit these data to an 


equation of the form 


Cc 
(22 eB eae 
oy 


n (6) (53-8) 


where n(6) = 5.56 - .0467 85 


and 85 is the Debye Theta 


3.8 at 86 bar 


3 
I 


nh = 3.6 at 120) bar. 
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Figure 31 shows a fit with the theory of Niklasson 
to crystal-I2.using the value of n from equation (5.3-8). 
Although this is a three parameter fit - two parameters 
to fit v and one to fit the knee - the relaxation t' 
is of significant importance in the region of the knee 
only. All the other parameters, except A and B in the 
adiabatic velocity, are experimentally determined. The 
parameter A is essentially normalized to unity by dividing 
by the isothermal velocity. Thus,we have, in reality, 
a one parameter fit. It is interesting to note that 
fee 1S wWithinuls+.OL Unity for the region & > 1.5k. 
I think this is sufficient evidence to state that the 
velocity above the knee is adiabatic. 

It is interesting to note that the (aP/aT) 
data of Jarvis et al (1968) shows a change in the power 
of the temperature change at about T/6 ~ .04. This 
change in power is reflected in their (c./c,, - 1) results 
giving rise to a small kink in the curves. The reason 


for this behavior is not known at present. 


Geetes Velocity Changevin ‘the Region 9 T= Te - Classical 
Theory 
The velocity change relative to the isothermal 


velocity in a classical model can be put in the form 


i 


c rae my 


one y 
7 
, 







— 


— 7 


-ai3 st dzgiw 312 s ewore If sivpl4 | ay, 
a® -« > F 5 Ts 5 a 
IDE 12 a 30 Suisv add patey SI isteyi> oF 


lire “ ~ ~ ~~ 4 rs —_ o “ 7 = " ; - i 

ry 2197emS1sq soTi? s ei aidz fpuodt lA 

a os pas pee al 
doa ‘t ony ~- Bgonn SAF II2 OF ono bas Vv. #23 oF 
. Cat 


I to aol; ont At sonetrogmi tnsocitinple Yo ek 










g 


A JG95xe ,@ie%emsA1sq Forljo eft LIA 
OS YALL SIGS TegGKS sis ,ytiscolev olisdseibs 


‘ Z Sitlemewan wo ff r #raheo.s » > 4 Sa —s 
Siis 7 (+45470S@!@9 2L A 19J9m 6 35! 
¢ ‘) 
: ay » bears Tce = wens 5 , 
2unT .yttoolsy Lsgarietitoel efs 


- 
sit oJ OF ssehive Jneioitive at eid? AdkAg ye 
mee 


(Isdeibs ¢i eend siz evods ysisofey 
is 


“> 

6. )\ 20s 35 spredo stTuys19eqmes of7 Yo 
Me 

pris 


13 mi 5estoaliay ad x9wog al s 






, 
a 
ee 
vir 

~b 







et ant .asviwo oid. ot Anid flame = oF seta 
| ra 


* > 
»3n9a91%qg 36 awond ton =2 xolverdsd naoads 10} 
U 
te eee 










o7 * 7 pokneh sits 
— 


1.000 


0.999 


0.998 


0.997 


V(T) + V(0) 


0.996 


0.995 


FIGURE 31 


Kiyo Ales 2 
86.4 BAR 
V (0)=669 m/s (428) 


cee CIE ORY 
rt THIS WORK 


1 2 3 
TEMPERATURE (KELVIN) 


Relative velocity change as a function of 
temperature for crystal I2 illustrating 
acaumacy of fit for 7 To 


118 





Resear Noe ep 








Bil 


Se mee ee 
; H | 


& 


rf 
wae 7 





$I sateva> 
AAG bOB 
(852) e\m 260-10) V 
YAO3IHT — 
AAOW 2IHT x 


PT 
= Was ss : 
— 


. 









where 


and f 


and 


To first order g 


or 


in 


A 
Q 
oS) Te 


N(T) is a~ l 


a Cc 
Vip (T) Kn? Vin CT) 
af? 
Nt eee cay 
f° + one 
is the frequency and 
cy" (T) 
£ = pn ES 
217 K uy 
S 
a=1l+ (~-1) ——+,—  . 
V Kn (T) op Vin (T) 


and hence 


119 


(Sao) 


efi 









[ 
(@-£.2) (D4 ———— a Soy args 
(7) pV yt 
S30 v0 
re (MD is 
ee 
Bbns yoneupss? oft et 2 
(t) “Wop 
Var “3 
ca) 1) 
—i——- -y+t=y + 
CT) pVqlT) yA Vv 


esoned bos i - 6 et (2) at » xSsbto texkd 


120 


Fon (I< Tor cee Stag the boundary scattering relaxation 


u B 
time. For our case Tp is, of course a constant, roughly 
5 





10 ~ sed. 
f in this region is 
aaa) 
(S| ele en Sem 
Vp 
Let 
Ty's 
RO gE ‘ 
u B 
thus 
Len 2 
(Tan 
MT) = 
2 ae 
Pr (5 ) 
in 
as 
imap 
3 v2 ; 
D 
Ps -22 
anus tom Tr 3-0; N(T) > (1/t,21) a elo for our case 
1 co, NIL) > 1 aS Tart, > On ke 


A rapid change should occur in N(T) about at 


Figure 32 shows crystal D4 fit to the above 


classical equations. Although one sees a small plateau 


OS1 









: - 
noijexsles pnitesiso= yxasbavod erly a? * 4? +57 >T 10% 


Yiripuoz ,inezenoo & gsexuOd Fo ,et a? e469 1y0 102 .emlis- 


2- 
ia 
2i noipez eid? al ~ 
s-" Ay g 
. i a? Re Ta 
av = 
ary? 
‘ > a’ 
a 
rs 
ee Oe 
( ) + Sy 
qi 


ss- 


8865 xU0 i012 Of « S (rSox\t) + (Tu .0 + f yxoR 


+ OF tig? Ge L+ (14 yee Tt \ 


36 duods (T)" at as 


V({T) +V(0} 


1.000 


0.999 


0.998 


0.997 


0.996 


0.995 


FIGURE 32 


CRYSTAL D4 
119.6 BAR 
V(0)=745 m/s (+35) 


—— THEORY (BHATIA) 
1 THIS WORK 





1 Z 3 
TEMPERATURE (KELVIN) 


Relative change in velocity compared to 
classical theory (Bhatia (1967)). 


“IS 





AG JAT2YaD 
AAG Self 
(2Et) em 2at=(O)¥ 


(AITAHS) YROSHT — 
ARROW 2IHT 2 










“ 
aq clip (LOG) 0). 20° t Ganeie a oeat? 













CRYSTAL D4 

119.6 BAR 

V(0)=745 m/s (+35) 
—— THEORY ( BHATIA) 
——THIS WORK 


0 1 | 2 3 
TEMPERATURE (KELVIN) 
FIGURE 33 Factor N(T) for crystal D4 compared to 


classical theory (Bhatia (1967)). 





Li22 


AG JATZYAaD 
GAB Sli 
(2E 2) vm @a0=(0)V 


{AiITAH4) YRO3HT ——— 


AROW C1HT—— © 













125 


at 1.7 K to 1.9 K, this is inadequate to explain the 
much larger plateau seen experimentally. 

A more severe test would be to compare the N(T) 
of the theory with the experimentally determined N, (T) 
from the above fit. This is shown in figure 33. 
While N(T) goes from 0 to l, N, (T) goesmirom -15i\ to 1. 
This is another sign of the weakness of the fit. This’ 
fit is typical of all crystals although the lowest N, (T) 
values vary from -36 to -8 depending on crystal and 
pressure. 

Thus, we must abandon the classical picture in 
explaining our results and go to a theory where N(T) 


can go negative. 


Bae «4 Velocity Change in the Isotropic Niklasson Picture 


There have been several theories dealing with 
sound velocity and its temperature change in the region 
Cat << 1 << Oty) notably (Gurevich and Efros, 1967) and 
the series by Niklasson, (1968, 1969, 1970,..1971). In 
the region of interest for our frequencies of excitation, 
these rieuries are essentially identical. 

Using the notation of Niklasson (1970) we have for 
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where s is the ratio of second sound to first sound. 


cs is the attenuation of second sound (with no coupling 


to lattice) 


PS = 5 ttt 57021 11 (Q)} 
u 

rs = 4 (+ + s?w*rt+(Q)} + s**t? (Q) 
u 


oo " ‘ ; 
where t -, tT , and t' are all relaxation times of the 


order of the total relaxation time Tt 


cel 
Ty Tp Ty 


For an isotropic crystal, rs becomes 
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Although the Niklasson theory seems to give a 
slightly better agreement for N(T), it is obviously 
still unsatisfactory. 

Figure 34 shows crystal D4 fitted to the iso- 
tropic Niklasson theory. A slightly better fit is 


possible if one varies Te and Tt however the fit is 


N! 
only slightly better. Figure 35 shows the factor 
N(T) with its experimental equivalent. The agreement 
here is qualitatively better. 

In general, Niklasson theory for the isotropic 
case is insufficient for our purposes. We can, however, 


parameterize the theory for a better fit. We will do 


this in the next section. 


gee -VeLocity Change with a Parameterized Niklasson 
theory 


Niklasson (1971) has attempted a parameterization 
of his theory for Argon, but our approach is somewhat 
different. 

As we need N(T) to become large and negative, 
varying 's' the ratio of second sound to first sound, 
will ene have a small effect. We must vary tT" (Q) and 
t"(Q). If we make 7T'(Q) 20 to 40 times larger than 


+" (Q) we can get much better fits. Figure 36 shows 
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erystal D5 with <° (0) = 23°* <"(0). Figure 37 shows 
the factor N(T). These fits are quite good. 

However, is there a physical reason for such a 
wide variation in t'(Q) and t"(Q)? Lets examine these 


relaxation times more closely. 
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gruneisen parameter and Vl ie Setne.maroscopic 


where y is the branch dependent (j) microscopic 
gruneisen parameter. Now the temperature dependence of 


y varies as much as 6%(Jarvis et al (1968)) and ina 


rather intricate fashion - thus there is probably some 
dependence on | q| of ee Asathe, c/agratio vio: He* 

is constant at most pressures, there is probably little 
orientation dependence in ap (Franck and Wanner (1970)). 
The vector nature of tT) (Q) will contribute to the varia- 


tion in this parameter. Whereas t"(Q) will be proportional 














eworle TE exupit .(O)"7 « ES = (O)'F doiw 20 Ieduyszo 
.boop stiup ots edFY eesdt . (T)H 103082 ont 

6 dove 102 nossez Isoleyig 8 sexsi? at , 19VEWOH. x 
seott onimaxe asad §(0)"7 baw (0)'r mi nOlssiusy ebiw | 


‘Yiseolo exom sami? nolssxslez © 


"1,000.0 % = = (0)*%r 


5ya8n byaw “O a 


syan” 7,9 Q,9 RS t= tort d 


‘ ” 
byan’ ~ byaa” ; 


das letayis 2iduo 10 god, s evel ew sa pool aa 
PP | 
2 (prt = (py Bey! 
r 


oigooesotsim (¢) tmebmeqeb donexd od at tty orertw 
oiqgoosoxsm eft al (p> bas tatemazag asatenre 

to sonebneqeb s1ts19qme3 eft wou . 18d ome 16q coarse 
S al bas ((8aeL) fs +s aseqesee 26 —s a6 sv yy 


omoz yidedoxg el sxedd eurid ~ m 






131 


to some average Tt - say rae) will be 


“isotropic! 
proportional to i}. Although TH is roughly isotropic 
ian hep Heres fs varies by about 50 depending on orient- 
ation. There will also be an effect due to anisotropy 
of x Ser of the order of the variation in velocity. -— 
about 203% 
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favors modes perpendicular to a, 
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u u 


These considerations seem to say that 117 (Q) 
will probably only be 5 timest"(Q) at these orientations 
at the most. 

Hence we must look elsewhere for this variation. 
If we examine the equation for the drift velocity 
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in the phonon drift velocity due to diffusion like processes, 
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lattice deformation. 

Now, recall that helium is a quantum crystal. 
Small changes in the local lattice will produce changes 
in the local lattice constants which in raen will affect 
the phonon spectrum. These changes in the phonon spectrum, 
in general, will mean quite long range effects in the 
lattice due to the consequences of the self-consistent 
nature of the phonon field. One can thus view Cyt! as 
a mean free path or coherence length which essentially 
gives the spatial dimensions of the effects of the local 
lattice deformation. It is not unreasonable to expect 
Cyt! to be larger in quantum crystals than in ordinary 
crystals, which is what is observed. 

We, therefore, believe that the large value of 
t' necessary to fit our data is direct evidence of the 


quantum nature of the solid. Notice that the variation 


of no other parameter could possibly fit our data. 


5.3 A Qualitative Description of the Attenuation 


‘During measurements, the attenuation was indirect- 
ly observed through the change in amplitude of the pulse 
being observed. The following gives a qualitative des- 
cription of our observations, and a comparison to the 
various theories. 

At Te T, we observed little attenuation (after 
seeing over 100 echoes), however when one neared Tor a 


rapid increase in attenuation occurred. This peak in the 
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attenuation was very broad and increased steadily to our 
lowest measured values. On one crystal (crystal D5) a 
small decrease in attenuation was observed at 0.8k, 
however, the decrease was very small. 

Figure 38 shows attenuation from the theories 
of Niklasson and Bhatia compared to our qualitative 
appraisal. Crystal D5 had about 70 visible echoes at 
3.6k, and 20 visible echoes at 0.8k - a change in 
attenuation of about 3.5. The theory of Niklasson 
appears to be in better agreement with our qualitative 
estimates, however, our measurements are much too 


primitive at this point to make any firm statements. 
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CHAPTER 6 
CONCLUSIONS 


4 


In single crystals of hcp He”, the velocity change 


in temperature has a sharp "knee" or anomaly at = ~ = 
This knee occurs where Ce ae l. It can be explained in 
terms of a parameterized Niklasson theory if we assume 
the second sound resonance predicted by Niklasson is 
much broader than the theory predicts. 

The high temperature region gives a good fit to 
the Bdasbaric velocity derived by Niklasson. The direc- 
tion of change of velocity in this region (T > tT) is 
decreasing for increasing temperature, lending support 
to the compressibilities calculated by Jarvis et al. 
(1968) over those calculated by Edwards et al (1965) 
calculated assuming a reduced equation of state. 

The attenuation shows qualitative agreement with 
Niklassson's parameterized theory, however, only rough 
qualitative data are available for analysis. 

Further work is necessary to determine the fine 
structure in the velocity plateau, as well as low tem- 
perature measurements to determine the value of the 
velocity in the zero sound limit. Measurements of the 
attenuation is necessary to determine which of several 


theories is more accurate. 
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APPENDICES 


Al Appendix for Derivation of Certain Thermodynamic 


Relations 
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A2 Niklasson Theory as Applied to Hexagonal Crystals 

Sound waves do not propagate perpendicular to 
the transducer in a hexagonal solid except at certain 
Symmetry points. Therefore, the phonon drift velocity 
in the crystal will in generral be at some angle 6(6), 
a function of orientation 6, to the externally applied 
field. 

Thus, if we excite a crystal at orientation 60, 
frequency & and wave vector Q, we will internally 
develop a phonon drift velocity at some angle 6+6 (98), 
with individual phonons of wave vectors and frequencies 
ae 


Assuming linear dispersion 


Q 


e, (8) (6) 
where &, (6) is the velocity (phase and group 
velocity are identical here) of sound mode teat: 


orientation 6, as c. and Q are colinear 
Be te AEM) s 


We will be performing averages over the internal 


variables 6', @¢' 
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A. (q,w) is the generalized spectral fit 


hiw) = , the bose occupation number. 
e 


If we assume harmonic phonons 
A; (aw) ~ AS (a1) = 4[6{wtw, (q) }+6{w-w, (q) 3] 


where te (a) = c,(8)q the harmonic phonon frequencies, 


We will be calculating averages over real quantities 


which depend only on q, es and x0!" 2 
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Section 2.7 gives an introduction to the theory 
by Niklasson wherein the integrals are calculated assuming 


an isotopic or cubic medium. 


~1 be our 


Now let 7(6') = eoagyanitiee — + oe 
B 


3 
Ty and Bp are the 


umpklapp, normal and boundary relaxation times, respec- 


total relaxation time where Tat 


tively. Only Te depends on orientation (assumed). 
As there is no dependence on Azimuthal orientation 
we can choose 6 = . and thus only x and 2z components 


will contribute. 


If we define 
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ee ey Gees eae 
cos 6 
and 
erat 4 
Sin 696 0 
e = 0 sin? Q' 
0 0 2 cos ou 


then we can write the solution to the Pieirls-Boltzmann 


equation (using the notation of section 2.7) as 
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for steady state 


iC . eRe 
= (po? 4. B+ G-B) 
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These calculations are in Q space. To convert to "r" 


Space one uses 


6 K 6 in q space is equivalent to 


ooo sin- 6 + — cos? abe 
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in real space. 

Now if we assume that coupling to lattice deform- 
ation adds only second order terms to the calculation of 
second sound, we can derive second sound by the following 


arguments. 
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where 1(Q,2) is the temperature gradient, 
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where Ch and Dy are the speed and attenuation of 


second sound respectively. 


In the normal region where second sound propagates 
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where 6 is the orientation of the crystal in real 
space. 


Define the generalized viscosity 
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Further, let the coupling between the temperature field 


and the lattice (which is proportional to O-W) be non- 
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We will now calculate the propagation of first 


sound with coupling to the local temperature deformation. 
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We must solve the simultaneous equations 


ek vege = <u> + = he<t> = - = + iW-QT(Q,2) 
[-iRC +0-K-Q]T(Q,2) = -pT2Q-W<u> 
where 
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sie 0 are the isothermal elastic constants 
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By Bay 
j = is the heat current 
W is roughly the thermal expansion. 


Solve the lower equation for T1(Q,®) 


T(Q,Q2) = eee 
<-inc +0-+K-G> 


and put this value into the upper equation 


youd wre Peta ecole 
p7<ue— 2 j-<ur+ 22 Feu> = - 2 3 SoTL (Owe (o-w) <u 
p p u Q [-iNC +6-K-O] 
let 
ee) {O°W W-O} + = = 
Fee ee ue pr pang haeo 
O-R:0 20 le) 
[-iNC_+0Q°K-Q]2 
then 


a~<w - = (W+¥) <> re oa (N+u)<u> =-= 9 


ele 


Thus X is the change in velocity and u the change 






enoitsups avosnsjiumie es svioae Jeum oW 


4 r ¢ 
, 7 af @ a { > 
be y [? oy > , <€o7> i— — <p> # 
= 4 G SO b~W I : 
on \ ow a : Ww A *O+ S2i- - 
Vv a 


imeiico Jsed eis et = + 
, fc Is) Ts | {3 i proy aft Ww 
~*~ . 
~ \ 1 > - — rf ‘ 
(i, t oJ ipa rawol oerlz 
. »f 7 ye - —_ - 
4 7 - -_ 2,0) T 


~ " 7 Ps 7 a tC , ba : 
<u>{(W-0) @(W*O) }Toi- < [ = {Le = | ¢ 
iti... C= > = <n == ¢<y>* = «<gs70 
[OrA+O+ ont—] Q _ : ba 





148 


in attenuation for first sound. 
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Projecting onto polarization vectors eJ (Q) and neglect- 


ing polization mixing 
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if the linear coefficient of compression is isotropic 
(Wanner (1970) ) 
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A3 A Method for Scaling Elastic Constants in hep He“ 

We have, at present, the elastic constants at 
four molar volumes from the work of Crépeau~ et al (1971) 
at 20.97 cc/mole, Greywall et al (1970) at 20.5 cc/mole 
ana rranek et al (L970) at 20232 and 19;28 cc/mole. 

These elastic constants were primarily determined 
from acoustic measurements. There was great variety in 
the method of determining the orientation, and some of 
the authors used various assumptions to calculate all 
the elastic contants. 

We would wish to calculate the elastic constants 
at any density, the following being a plausible way of 
going this. 

Pirst, Franck et al (1970) proved that as.a 
result of the c/a ratio remaining constant in hcp He4 


to very high densities 


€ ne =C + € ° (A3-1) 


Using this equation one can simplify the expres- 


sion for the compressibility to 
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The compressibility has been measured from (SF) 


by Warvis 6 61) (196e) trom l/.77 to, 20. /2 cco/mole. 
Recalling that the Debye Theta has also been 


measured by Ahlers (1970) at a wide range of molar 
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volume, one can use these data to scale the elastic 
constants. 

First we must modify the data of Crepeau et al 
(1970) to fit condition (A3-1). This can be done within 
their experimental error. 


Then at all four molar volumes 
a— = 1.34.01 . (A3-3) 


It should also be noted that in this form systematic 
errors in the measurements (such as distance calibration 
error, or time scale calibration error) are negated. 

Also C33 and Cia are probably the simpliest elastic con- 
Stants to find as they are simply related to the Long 
itudinal velocity along the 'C' axis and in the basal 
plane respectively (see section 2.3). 


Recall, also, that 
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or the transverse velocities are directly proportional 
to the Debye Theta. 


Further let 
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We will further assume that the ratios of the 
transverse velocities o remain constant at all densities 
2 
Caa 


We can finally write 
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From experiment o is in the range 1 to 2. If one chooses 
a value of o and then calculates A to give the correct 
On! One can get reasonable fits. The final fit is 
Surprisingly insensitive to the size of oa. 

In extrapolations of elastic constants for com- 


parison with theories, I have used 
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and 


¥ cg 270 +83 (V-14. 208) 


Gow =.168 56 {sa—a08 


This gives Kn to 1% and or tor 

This scaling gives results within experimental 
error, except for the data of Crepeau et al (1970) 
which appear about 2% low. Also, this scaling gives 


results consistent with the velocities I have measured | 


at 120: bar-~and. 86+ bar. 
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A4 Computer Programs 

This appendix gives some of the computer programs 
used during the completion of this work. 

TEMP is the function that converts resistance to 
temperature. COTEMP and COXTEMP are the coefficients 
used for the main region of fit and for extrapolation 
respectively. LL and LU define the main region of fit. 


PO fits a number X to a polynomial with coefficients 


VPLOT calculates and plots the velocity versus 
temperature from the raw experimental data. 

VSS is used to give the printout in the next appendix. 

It fits the velocities from VPLOT to the theories of 
Bhatia or Niklasson. 

CODEX arranges the raw data. SCALE calculates the 
elastic constants Ae a given pressure. THETAV calculates 
the Debye theta for a given volume, while KTVM calculates 
the compressibility. DEBH calculates the Debye theta from 
the elastic constants for hexagonal crystals. BHAT is the 
theory of Bhatia , while VEL gives the theory of Niklasson. 
SDFIT, FIT, WFITL, and LINEAR are functions that calculate 
the least square coefficients. REL calculates the relaxa- 
tion times. UMKLAPP calculates Umklapp relaxation times. 
CV calculates the heat capacity. ORVL calculates the 
orientation. PVDEG produces a polynomial of fractional 
or negative powers. 


HEXA gives the velocities and related quantities for 
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hexagonal crystals. DELL gives the deviation from a pure 
mode. DELTAL, DELTAT1, DELTAT2 give the deviation from 
wave normal. VL, VTl, and VT2 give the velocities. PH is 
a subroutine. 
VEL2,INSET, and INFUN are functions to calculate 
the orientation dependence of the velocity of second sound. 
Finally, DEBC, VCUB, CUDEF, and RCR are used to 
calculate the Debye theta of cubic materials. RCR is a 
functions that gives the three real solutions of a cubic 


equation if they exist. 
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A5 Experimental Data 


This appendix gives the output of function VSS. It 
is a compilation of the experimental data on twelve crystals 
and the fits to these crystals fron the theory of Niklasson. 
The headings give the temperature, normalized velocity, 
attenuation, and N(T) factor respectively for the theoreti- 
cal fit. The next page gives the experimental data for 
the same crystal. The headings give the temperature, 
absolute velocity, normalized velocity, and N(T) factor 


respectively for the experimental data. 
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TSOTHERHAL VELOCITY = A+F 
A27TU209 + /FL 8. 38 M/S 


CRYSTA D403 


DERYE @=42.76 


SECOND SOUND=206.1 MIS 
SLANDARD DEVIATION FOR FIS@ES 
RELAXATION FACTORS 1xoH 
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0.993848 
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1.000 
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1.009 
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1,000 
1.000 
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EXPERIMENTAL DATA FOR CRYSTAL 03 
vi riev[OJ 
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ISOTHERMAL VELOCITY 
A=750.8 +/- 32.84 M/S 


CRYSTAL A8 


DEBYFE @=42,25 


SECOND SOUND=203.4 M/S 
SLANDARD DEVIATION FOR FIT IS 


3120.8 BAR 


AtBxi«l 


330/08/71 


B= 0,11139 


MOLAR VOLUMFE=17.439 CCO/M 
FREQUENCY=S5E6 HZ 
0.08358 


RELAXATION FACTORS 1.5xTN,3.7xTU,6xTAUL 


ORIFNTATION IS BETWEEN 40 AND 30 


THEORETICAL CALCULATIONS (NIXLASSON) 
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0,99722 
0.939685 
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0.99308 
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0.575 
0.693 
0.825 
0.967 
3. BHO 
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1.432 
1.477 
1.480 
1.440 
1.356 
1 236 
1,088 
0,930 
0.785 
9.662 
04902 
0.478 
0.409 
0.392 
0,304 
0,265 
0.232 
0,204 
0.181 
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7,884 
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0.744 
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0.967 
0.983 
0,991 
07925 
0.998 
09.999 
04393 
1,000 
1.000 
1,000 
1.000 
1,000 
1.000 
1,000 
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RPXPFERIMENTAL DATA FOR CRYSTAL, AB 
V(UrTJ+evl[o] 


TEMP 
3.530 
3.464 
3.386 
3.233 
3.243 
3.183 
Os 
Aeo82 
2696 
2,804 
2.708 
2501 
a. 318 
2.210 
25, 096 
i772 74 
19,609 
1.413 
1 262 
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vo mJ 
toss Cee ty 
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753.48 
WSS7S7 G 
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754.73 
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CRYSTAL C2 


ISOTHERMAL 


DFRYF @=42 


SFPCOND SOUND=202.5 M/S 


tio eo BAR 


VELOCITY = 
A=738 +/- 43.18 M/S 


. 09 


At+BxT un 


B= 0.1185 
MOLAR VOLUMF=17.52 CC/"™ 
PREQUFNCY=5F6 HZ 


STANDARD DFVIATION FOR FIT IS 
RELAXATION PACTORS 1.5x7NV,5.4x7U 
ORIENTATION IS BETWEEN 45 AND 30 


-O1 
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THEORETICAL CALCULATIONS (BHATIA) 


TEMP 
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0.953942 
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0.99776 
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0.99735 
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0.955 
0.975 
0.985 
0,991 
0.39395 
0.997 
0.998 
0.999 
0.993 
0.999 


N=3,592 
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EXPERIMENTAL DATA POR CRIOTAL C2 


TEMP VOTs Vi TI+VC0J NIVPAC 
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3.393 733°329 0.39349 4, 002 
3,309 7333 $60 0.993393 9.990 
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ISOTHERMAL VELOCITY = A+BxT«N 
A=758.2 +/- 22.75 M/S 


CRYSTAL D3 


DEBYF @=42,.,0b 


SECOND SOUND=202.3 M/S 
SYANDARD DEVIATION FOR FIT IS 
RELAXATION FACTORS 1.5x7N,6.8xTU,4SxTAY 


ORIENT, 
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B=70.1025 


N=3.594 


MOLAR VOLUMR=17.52 CC/M 


40 AND 55 OR 65 AND 30 
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EXPERIMENTAL DATA FOR CRYSTAL D3 
VO TitV[OJ 


LEMP 
3.6193 
3.568 
3.484 
3.414 
3.361 
oeod 
8,151 
3.060 
2.951 
2.846 
2.840 
2.692 
2.549 
2.308 
1.608 
fexvas 


gral 
753.94 
794409 
TS 4, Ov 
T5456? 
754.96 
758729 
SIG DF 
755.85 
756.20 
4564.55 
756263 
TSOFES 
756.63 
750.63 
756303 
750,63 


0.99440 
0.99459 
0.99483 
0.99536 
0.99574 
0.99618 
0.99655 
0.99691 
0.99737 
0.99783 
0.99794 
0.99794 
0.99794 
0.99794 
0.99794 
0.99794 


NIPPAC 
1.004 
0.993 
0.962 
0.986 
1.006 
1.002 
0.989 
0.990 
Olt 
17.054 
NRO WAY | 
0.948 
0.730 
0.410 

92.3973 

342.103 
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ISOTHERHAL VELOCITY = AtRxT«N 
B= 


CRYSTAL Dy 


A=754.1 +/- 4.223 
DEBYF @=42,.1 
SECOND SOUND=202.6 M/S 


STANDARD DEVIATION FOR FIT rs 


3103. 6 


MIS 


BAR 


t2Sh Ia 7/74 


0.1161 
MOLAR VOLUMP=17.52 Ccys™ 


Vege 591 


PREQUENCY=SF6 #Z 


9.04922 


RELAXATION FACTORS 1.5xTN,2xTU, 30x TAUL 


ORTENTATIO”N IS BETWEEN 4O AND 90 


THEORETICAL CALCULATIONS(NIKLASSON) 


TEMP 
0.3 


e 
MON OMP OD OYHUME OHH OU NOM FONE Ss 6S SOG Se 


virji+vCoJ ATTEN 
0.99994 0.261 
0.99984 0.379 
0.99966 0.502 
0.99937 0.653 
0.99300 0,859 
0.99863 16,137 
0.99837 1.474 
0.99830 12324 
0.99838 25.120 
0.99856 2532.9 
0.99876 2.426 
0.99890 on 04 
0.99893 226A. 
0.99884 2. 022 
0.99867 155704 
0.99853 1.329 
0.99852 0.973 
0.99861 0.709 
0.99864 0.539 
0.99858 0.429 
0.99844 0.350 
0.99824 0,290 
0.99799 0.244 
0.93770 0,207 
0.99737 Oh tae, 
0.99701 Oni5e 
0.99661 0.132 
0.99617 Cr he 
0.99568 0.102 
0.99515 0.091 
0.99458 0.081 
0.99396 0.073 
0.99328 0,066 
0.99255 0,060 
0.99176 0.955 


NIVFAC 


“50,654 
“50,369 


43.462 


Page ou 
“42,820 
“35.956 
et ee 
“49.199 
128.396 


“72.350 
“4.098 
725 1:78 
7415 193 
70.746 
“0,543 
08 303 
0,084 
0.476 
0.737 
0.874 
0,940 
0.974 
0.985 
0,992 
0.996 
0.998 


07. 999) 


0,993 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
1.000 
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FPXPERIMENTAL DATA FOR OR YS TAL: Ph 
VE 7IJ+VL0] 


TEMP 
3.552 
3.623 
34567 
3.495 
3.419 
3.358 
3262 
3212 
3.070 
2.988 
2e 991 
2721 
22380 
2.406 
2.164 
13571 
1.360 
0.763 


ViTJ 
TROD 
748.33 
748,61 
749.03 
743.45 
749,79 
750,26 
750.44 
750,33 
78420029 
TSiTF74 
752207 
752.43 
TS2372 
753.12 
753.16 
753.18 
753 922 


0.99205 
0.99230 
0.99272 
0.99329 
0,99384 
0.99429 
0.99492 
0.99515 
OF 99584 
0.99628 
0.99688 
0.99731 
0.99780 
0.99818 
0.99871 
0.99876 
0.99879 
9.99883 


NIVPAC 
0.987 
0.998 
0.990 
09,996 
£004 
O24 
1% 622 
4.012 
0.995 
1.016 
0,997 
i004 
7. 006 
0.951 
6.910 

S25129 

wo .987 

“30,490 
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CRYSTAL D5 


$113.5 BAP 


ISOTHERMAL VELOCITY = A+BxTuN 


A=757.9 +/- 13.49 M/S 


DFEBYE @=42,1 


SECOND SOUND=202.6 M/S 


so0/19/71 


B= 0.1147 


N33 ,591 


MOLAR VOLUMF=217,52 CC/M 


SLANDARD DEVIATION FOR FIT Is 
RELAXATION FACTORS 1.Sx2NS165xTU, 33x TAU 
ORTFNT. IS BTWN. 40 AND 55 OR 65 AND 90 


THEORETICAL CALCULATIONS (NIKLASSON) 


TEMP 
O73 


o- OF 6M OF om eh) oh off « 
YOON FWNHRPOUWDYIA UNE 


PRRPRPRR OOOO COO 


NOWNFWONRP OW DIMNHVE WHHR OO Oo 


WOWWHWWWWHNNHKYNH HNN PDB Epp 


VCOTItVL OJ ATTEN 
0.99994 0,260 
0.93983 0.378 
0.99962 0,502 
0.99937 0.657 
0.99890 0.872 
0.99850 1.166 
0.99822 1.524 
0.99815 1,890 
0.99826 2,495 
0.99849 2,009 
0.99872 2.444 
0.99887 2.357 
0.99889 2.146 
0.99878 1.838 
0.99863 1.460 
0.99858 4.066 
0.99807 0,747 
0,29877 0,941 
0.99877 0.415 
0.99868 0.330 
0.99853 O39259 
0.93832 de2eo 
0.99807 0,186 
0.99779 Os Lo 
0.99747 0.134 
0.99712 0,115 
0.99673 0.100 
0.99631 0.087 
9.99584 0’. 077 
0.99533 0.068 
0.99478 0.061 
0.99418 0,055 
0,99353 0.050 
0.99282 0,045 
0.99206 0,041 


NIKPFAC 
“55,806 
“55.492 
754.495 
“52.033 
"47.190 
“39.630 
"30.308 
"21 063 
"13.302 

“7,853 

“4.300 

72 DOO 

Lt1302 

“0.915 

“0,653 

0. 2'59 

O 237 
0.614 
0.820 
W207 
O.954 
0.981 
0.991 
0.995 
0.997 
0.999 
0,999 
1.009 
1,000 
1.000 
4 0010 
1,000 
1,000 
1,900 
4,000 


PREQUFNCY=SE6 AZ 
0.0421 
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EXPERIMENTAL DATA FOR CRYSLAL DS 
V(TJ+VC0) 


TEMP 
3.638 
3.600 
3.553 
3.487 
3.423 
37350 
3.278 
3.194 
3.106 
3.010 
2895 
20749 
25607 
2.4539 
2.403 
2.230 
2.109 
Le SR7 
1.370 
0.662 


Vit J 
752,20 
732, 39 
752.66 
733.90 
753,33 
753.68 
754,04 
754,38 
TORTI 
755.04 
755930 
7553.82 
756.18 
756.46 
756.57 
756.86 
756.89 
756,91 
756,90 
756.90 


0.99253 
0.99279 


-0.99314 


0.99358 
0.99402 
0.99449 
0.99497 
0.99540 
0.99584 
0.99628 
0.99675 
0.99729 
0.99778 
0.99816 
0.99830 
0.99868 
0.99871 
0.99874 
0.99873 
0.99873 


NIXFAC 
0.939 
0.995 
0.998 
0.996 
0.998 
1,001 
1.011 
1.007 
1700G 
1.004 
0.937 
0.994 
a. .008 
0,978 
0.976 
0.956 
0.785 
0.048 

63.215 

"67.525 
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ISOTHERMAL VELOCITY 
A=742.1 +/- 13,28 “7/3 


CRYSTAL F2 


DFBYr @=41,.88 


SPCOMD SOUND=201.4 M/S 
SLANDARD DEVIATION FOR Frr rg 


$117..8 RAR 


A+BxTxN 


20 fA 2/74 


R= 0.1162 


RELAXATION FACTORS 1.2% DN 48 
ORIFNTATION IS BFETWPREN 4O AN 


MOLAR VOLUMP=17,55 CCA 

FREQUPNCY=5E6 HZ 
0.04584 
OxTU,27xTAUL 


THEORETICAL CALCULATIONS (NIKLASSON) 


TEMP 
0.3 


PrPPRPPRPPRP PEP PRPOGCOCCO 
e e e e e e e e e e es e e e ° e e 
ONFWNKrFOWOYCHOE 


CEP aeen tO Ole” Cel 6) ee Med eo 6.) 
PNMNEFWNHFOANW DUM EWHHROWOY 


WWOWHWWWWWWNHNDHNYN HO NNNDD\h 


VO TJeVL0 ol 


0.99995 
0.993986 
0.99968 
0.99942 
0.99908 
0.99874 
0.99851 
0.99846 
0.99855 
0.99874 
0.99892 
0.99903 
0.99902 
0.99891 
0.99876 
0.99871 
0.99877 
0.99882 
0.99877 
0.99865 
0.99847 
0.99825 
0.99798 
0.99768 
0.99735 
0.99697 
0.99650 
0.99611 
0.99562 
0.99508 
0.99450 
0.99386 
0.99317 
0.99243 
0.99162 
0.99076 


ATTEN 
0.270 
Oe oon 
0.917 
0.668 
0.871 
1.142 
1.466 
Le. 7 92 
2.001 
2,,228 
2.269 
7 es Gl 
ae oe 
1,695 
1,348 
0.994 
Ore ed 
0.528 
0.419 
0.329 
QO. 269 
0.223 
0.186 
O,1o7 
9.134 
0,115 
0,100 
0.088 
0.077 
0,969 
0.061 
0.055 
0,950 
0.046 
0,042 
0.038 


NIKFAC 
“45,499 
SSE Oat 
TH. 394 
the as oO 
“38,259 
731.959 
"24 955 
"10,684 
“10.460 

Beene 

5.30,7 a 

oie 96 

0,848 
“0,536 
“0.319 

0.019 

0,420 

O24. t 

0,867 

0.9393 

Out 

0.986 

0.993 

0.9956 

0.998 

0,999 

0.999 

1.000 

1.0006 

1.000 

1% 06.0 

1.000 

1.000 

1,000 

1.000 

1,000 
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EXPERIMENTAL DAT 


LEMP 
3.726 
3.692 
3.653 
3.594 
3.492 
3.2433 
3.9323 
< erate 
3.174 
3.066 
2: DF 
2.838 
2.696 
265 GS 
2.444 
Bienes dal 
$997.5 
Tal tad 
260-15 
1.368 
0.641 


VOT 
Cis ergy? 
735.93 
736.18 
736.54 
a Wy ee | 
737.44 
737.96 
738.30 
738.63 
739, 20:7 
739.38 
739...79 
740,17 
740.43 
740.70 
740,99 
7356 2.2 
7 e502 
741.33 
tS 
744.33 


A POR CRYSTAL F2 
VOTJevCo] 


0.99195 
0.99103 
0599197 
0.99246 
0.99323 
0.99368 
0.99438 
0.99482 
Oe a5 2:7 
0.99587 
0.99629 
0.99684 
0.99734 
0.99770 
0.99807 
0.99845 
0.99877 
0.99886 


0.99891 


0.99891 
0.998914 


NIVFAC 
0.995 
0.993 
9.995 
0.999 
1.000 
1.005 
1,009 
1.009 
1.0093 
1-016 
1.009 
1.005 
0,992 
0.985 
0.959 
0.925 
0.609 
0,141 

PME TYTE 


_ 2,340 


63,089 
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ISO?THFRMAL VELOCITY 
A=718.4 +/- 139.11 M/S 


CRYSTAL F2 


DFERBYF @=42,03 


$119.0 BAR 


A+BxTuil 


2301/02/72 


B= 0.1013 


MOLAR VOLUMPF= 
SECOND SOUND=202,2 M/S 
SLANDARD DEVIATION FOR 
RELAXATION FACTOR 
ORIPNTATION IS Br 


FIL IFS 


gM A ad eo 
PREQUFENCY=5E6 HZ 
0.04307 
S 1.5x?V,0.3«xT7U,17xTAU1 
PEN 50 AND 70 


THEORETICAL CALCULATIONS(NIKLASSON) 


TEMP 
O23 


on) oy.e 
one 


ee 
on) CC: 


VNP RP PRP PPP HE GO OC OC oO 
e 


WWWWWWW WHY DY HNN PL 
e 


VE7I+VL 0] 


0293937 
0,99994 
0.994980 
0.99464 
0.99942 
0.99921 
0.99909 
0.99942 
0.99926 
0.99944 
0.99951 
0.99949 
0.99946 
0.94948 
0.99949 
0.99946 
0.99938 
0.99927 
0.93913 
0.39898 
0.99880 
0.99861 
0.99838 
0.99814 
0.99786 
0.997556 
07939723 
0.99687 
0.99647 
0.99604 
0.99557 
0.99506 
0.99450 
0.993990 
0, 39925 


ATTEN 
0.275 
0,397 
0.520 
0.558 
0.826 
1,034 
1.261 
1.443 
1.4398 
ao od 
Sedo? 
0.896 
0.631 
0.422 
0.292 
0,213 
0.160 
7223 
0.096 
0.076 
0.061 
0.050 
0.041 
0.035 
05023 
0.025 
0.022 
0.019 
0.016 
0.9015 
0.013 
0.012 
0,011 
0,020 
0.009 


NIXFAC 
e26.308 
L2s 268 
TOPE 4E 
P26. 361 
729822 


7.19), 78s 


ise Sos 
“9.085 
“4.567 
T4810 
“0,589 
“0,135 
0.182 
02654 
OX 812 
0.929 
0.973 
0.990 
0.996 
0,998 
0.999 
1.000 
4°3000 
12000 
12000 
12000 
sou 
000 
1.900 
TOO 
1.000 
1.000 
Pe o00 
1.000 
1.000 
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FXPPRIMENTAL DATA POR CRYSTAL F2 


TH Mp YE7] VOTIJ+VCOJ PIKPAC 
3.637 713.83 0.99303 0.996 
3.621 7 1587.93 0.99377 A¥50:010 
34561 714,13 0.99405 0.988 
3.558 714.33 0.99433 Imi 2 2 
3.423 FAM i7 2 0.99487 0.996 
3.336 715.09 0.99539 1.006 
3.225 715.45 0.99589 0.994 
3.166 TAS O72 0.99627 Pots 
3.022 716.06 0.99674 0.989 
2.945 76a D7 0.99703 0.989 
2.687 716.88 0.99788 0.994 
2.539 P47 0.99829 0.999 
2, 273 T4179, 53 0.99879 0.958 
2.054 717.84 0.99922 1.014 
1.868 718.04 0.99949 1.066 
1.418 718.04 0.99949 “0,063 
0.640 718.04 0.999u9 "29.706 
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ISUTHERMAL VELOCITY 
A=732.4 +/- 11.72 M/S 


CRYSTAL PF3 


DEBYE e=42.06 


SECOND SOUND=202.3 M/S 
SLANDARD DEVIATION FOR FIT te 
RELAXATION FACTORS 1.5xTN,0.7xTU,23xTAU1 
ORIENTATION IS BETWEEN 45 AND 85 


2113.2 BAR 


At+BxTuN 


313/02/72 


B= 0.1046 


MOLAR VOLUMF=17.52 CC/M 
FREQUENCY=5E6 HZ 
0.05129 


THEORETICAL CALCULATIONS( NIKLASS(OM) 


TEMP 
0.9 


o ¢$e¢éeé@6mhCcithC<iC HC CO ]UCOS 
womMaMnn £ 


eyie ie. on. 6 
aoFWNFrO 


ee 
b> me?) 


p> rPRPPRP PPP RRR OOO OC Oo 
® 
oe) 


e e e e e e e e e °° ° e 


WWWWHWWNHNHHYDHYNH PDP NHNNA 
e 
ONFWNRP OW DIMA EWNHOYW 


VI TI+VL0] 


0.99996 
0.99988 
0.99973 
0.99954 
0.99923 
0.99894 
0.99876 
0.99873 
0.93886 
0.99906 
0.99923 
0.99928 
0.99923 
0.99915 
0.99916 
0.99923 
0.99924 
0299976 
0.99907 
0.99893 
0.99876 
0.99855 
0.99833 
0.99807 
0.99779 
0.99748 
0.99714 
0.99677 
0.99636 
0.99591 
0.99542 
0.99489 
9.99432 
0.99370 


ATTEN 
0.269 
0.390 
0.514 
0.658 
0.846 
1.090 
1.374 
1.646 
1.837 
ans 4 
1.796 
1.582 
1,295 
0.976 
0.683 
0.476 
0.349 
0.268 
Oe22:0 
0,168 
0.136 
Onair 
0.093 
0,078 
0,066 
0.057 
0,049 
0,043 
0.038 
0.033 
0.030 
ORO 27 
0,024 
6,022 


_NIKFAC 

38.633 
Poe vei6 
7.37 a? Le 
“35.983 
“32,578 
2272230 
“20,567 
"43° 2820 


“8,169 


“4 219 
~1,945 
0,884 
70,455 
“0.165 
0,229 
0,506 
0.828 
0.928 
0.970 
0.987 
0.994 
0.997 
0.999 
9.999 
12000 
1.000 
1°30 010 
1.090 
1.000 
1.000 
1.000 
1.000 
1.000 
P0059 
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EX PERIMENTAL DATA FOR CRYSTAL F3 
VOTlJ+vC0l 


TEMP 
3.594 
PSs 
3.431 
34231 
3.114 
3.002 
2.940 
27a 4 
22693 
2.470 
2.308 
2.473 
1.948 
1.807 
1.360 
1,349 
LL. 
0.7156 
0.655 


VOT 
727,85 


728,10 
728.55 


729,28 
729.869 
730.05 
730.20 
7390.48 
730,85 
783405 
731.43 
731,78 
731.87 
731.87 
731,87 
731.87 
731.87 
234687 
731.87 


0.99378 
0.99413 
0.99473 
0.99573 
0.99629 
0.99678 
0.99699 
0.99736 
0.99787 
0.99829 
0.99866 
0.99914 
0.99926 
0.99926 
0.99926 
0.99926 
0.99926 
0.99926 
0.99926 


PIVEAC 
1,005 
1.905 
1.001 
0.994 
0,999 
1.004 
0,999 


0.943 


1,020 
0.945 
0.946 
i445 
0.934 
0.678 
S40 245 
71,333 
“1.650 
“284324 
“39.770 
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CRYSTAL G1 3 86.2 BAR 315/04/72 


ISATHRERMAL VELOCITY = A+BxTuN 


A=754 +/- 13.04 M/S B= 0,1625 N=3.807 
DPBYF @=37.48 MOLAR VOLUMF=18.34 cCy/™ 
SECOND SOUND=178.9 M/S PREQUFNCY=5E6 RZ 


STANDARD DEVIATION POR FIT IS 0.0142 
RELAXATION FACTORS 1.5xTN,1xTU,4xPAU1 
ORIENTATION IS BRTWRER 10 AND 30 


THEORETICAL CALCULATIONS (NIKLASSO?N) 


TEMP VLTJ+#VL 04 ATTEN NIKFAC 
0.3 0.99999 0.384 6,001 
0.4 0.939997 0.540 “5.932 
0.5 0.99992 0,692 arts 
0.6 0.99986 0.850 S199 
0.7 0.99978 1.09 “4.254 
0.8 0.99971 1.196 "2,981 
0.9 9.99969 1.358 “4.653 
Deo 0.99970 4.870 “0.554 
tint 0.99973 1.497 0.184 
133 0.99975 1.424 0.583 
q53 0.99972 13266 0.750 
14 0.99967 1.057 0,810 
4,5 0.99360 0.835 0.855 
£26 0.99954 0.638 0.916 
Ty 0.99947 0.486 0.963 
1.8 0.99938 0.374 0.985 
1.4 0.99925 0.291 0.994 
Papas, 0.99909 0.223 0.998 
Ziel 0.99891 0.183 0,999 
es 0.99870 0,148 1,000 
edo 0.99845 Gatos 1.000 
2.4 0.99817 0.100 1.000 
BAS 0.99786 0,084 1.000 
226 0.99751 0.072 1.000 
fay 0.99711 0.062 1,000 
226 0.99667 0.053 1.000 
2.9 0.39617 0.047 12000 
a0 0.99562 0,041 t. 000 
Sat 0.99501 0007 1.000 
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FXPFRIMENTAL DATA FOR CRYSTAL (1 
vVCTJ+vCO) 


TEMP 
3.047 
3,012 
2.974 
2.914 
2.840 
2.753 
2.649 
2.543 
2.420 
ae aus 
2,092 
1237-0 
AV 726 
L677 
1.523 
13339 
1.150 
0.657 


VOT 
750.45 
750.63 
750,80 
Bal 02 
Tae a eo 
752.66 
75.2%99 
752,29 
YEP | 
752.93 
753.18 
753.41 
753,53 
753.62 
753.568 
753,74 
753.81 
753.81 


0.99531 
0.99554 
0.99577 
0,99613 
0.99649 
0.99690 
0.99734 
0.99774 
0.99812 
0.99869 
0.99893 
0.99922 
0.99939 
0.99950 
0.999593 
0.99967 
0.99975 
0.93975 


NIVFAC 
0.996 
0,999 
1.003 
POOL 
Ii .012 
1.004 
1.005 
17,004 
0,999 
1.000 
0.989 
6.9074 
0.967 
0.874 
0.706 
0.580 

“6,649 
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ISOTHERMAL VELOCITY 
A=674.7 +/- 6,207 "“/S 


CRYSTAL H2 


DEBYE @=37.54 


SECOND SOUND=179.2 M/S 
STANDARD DEVIATION FOR FIT IS 


36,6 BAR 


A¢ixTaN 


3904/05/72 


3200-1557 


MOLAR VOLUMF=18,32 CC/M 
FREQUENCY=5F6 HZ 
0.06896 


RELAXATION FACTORS 1,.5x2N,1xTU,1xTAU1 


ORIFNTATION IS BETWREEN 40 AND 30 


THEORETICAL CALCULATIONS (NIKLASSOM) 


LEMP 
0.3 


WNHNHHNYNNDNHNNN PRP PBR PRP PREP R COC COO 
PWDNONFWNKrRPOWDNIMDHEWNHYKrPOWDONOM Ee 


VUTJ+vLo] 


1.00000 
0€929999 
0.99998 
0.99995 
0.99993 
0.99990 
0.99988 
0.99986 
0.99985 
0.99982 
0.99977 
0.993971 
06399962 
0.99951 
0.99939 
0.99924 
0.99907 
0.99887 
0.99863 
0.99837 
0.99806 
0.99771 
0.99732 
0.99688 
0.99639 
0.99584 
0.99523 
0.99455 


ATTEN 
0.427 
0.599 
0.704 
0.928 
£.093 
1.253 
2.3933 
1.485 
1.503 
1.432 
1.282 
1.086 
0,880 
0.6393 
0.539 
0.418 
0,326 
0,257 
0,205 
0.166 
0.136 
0.113 
0.095 
0,081 
0.069 
0,060 
0.053 
0,047 


NIKFAC 
0.855 
“0,836 
m0 ST 
“0,636 
“0.381 
"0,035 
0.328 
0.628 
0.8293 
0,938 
0.982 
9.993 
0.993 
0.994 
0.997 
0.999 
0,999 
1.000 
1.000 
1.000 
1.000 
1.000 
1,008 
4006 
1.000 
1.000 
1.000 
1.000 
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EXPERIMENTAL DATA FOR CRYSTAL H2 
vO 7J+VC0] 


TEMP 
2.969 
2.965 
2.940 
2.894 
2.847 
2.808 
Zune tay 
2.749 
2.686 
21562 
2.551 
2.526 
2.486 
2.395 
mas S 
2.252 
2214 
2..%B 4 
21-246 
2.064 
1.986 
1.941 
1.3909 
1.864 
1.995 
1.961 
1.930 
1.850 
1.586 
1.704 
1,545 
1, 41 
19.225 
e187 5 
0.647 


vo nd 
Gee ig 
671.193 
G71, 129 
671.46 
671,55 
671,82 
671.96 
672.05 
672,26 
6072.65 
672,68 
672,73 
672,80 
673.04 
673,27 
673.39 
673.48 
673.54 
673.66 
673.74 
673, 87 
673.92 
673,36 
674,02 
674,10 
674.12 
674.15 
674,21 
074,24 
074,29 
674.42 
674.45 
074,50 
674,50 
674,59 


0.99482 
0.99487 
0.99500 
0,99526 
0.99554 
0.99581 
TuaI3GOr 
0.99614 
0.99645 
0.99703 
0.99707 
0.99714 
0.99733 
0.99761 
0.997395 
0.99813 
0.99826 
0.99835 
0.99853 
0.99864 
0.99883 


- 0.998914 


0:..9989 7 
0.99905 
0,99917 
0.99920 
0.199925 
0.99935 
0.99939 
0.99947 
0.99966 
0.99370 
0..19:9:9:7.8 
0.99978 
0.99978 


NIMPAC 
1.006 
1.008 
1.004 
0.999 
0.998 
1. 004 
1.004 
1.002 
0.998 
0.994 
0.993 
0.986 
0.990 
0.971 
0.983 
0.975 
0.979 
0.976 
0.974 
0.967 
0.970 
0.956 
0.951 
0.950 
1.143 
1.130 
1.128 
4.422 
0.837 
4.074 
1.068 
4.004 
0.846 

bi. 7 25 
6 292 
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ISOTHERMAL VELOCITY 


CRYSTAL I2 


DFRBYE @=37,51 


SECOND SOUND=173.1 M/s 
OLTANDARD DEVIATION FOR FIT IS 


86.4% BAR 


A+BxTuN 


32 405/ 72 


B= 00,4612 


MOLAR VOLUME=18,33 ccys™ 
FREQUENCY=SE6 RZ 
0,05547 


RELAXATION FACTORS 1.5xTN,0.2xTU,8xTAU1 


ORIENTATION IS BRITWEEN 45 AND 90 


THEORETICAL CALCULATIONS (NIKLASSON) 


TEMP 
0.3 


ee e e e ° e e e ° e e e e e s e 


DNHONNNYNNNAHNNNHP RFP RPP RRR HBHFRP OCC OOO 
e 


vCrjevloJ 


0.99998 
0.399993 
0,99985 
0.99972 
0.99957 
0.99947 
0.99950 
0.99902 
0.99968 
0.99966 
0.99962 
0,99960 
0.99954 
0.99944 
0.999314 
0.99914 
0.99895 
0.99872 
0.99846 
0.99816 
0.599782 
0.99743 
0.99699 
0.99649 
0.99594 
0.99532 
0.99463 
0.99387 


ATTEN 
0.433 
Of 670 
0.787 
0,979 
13192 
14.393 
1.510 
1.421 
1,158 
0,844 
0,565 
0,369 
On. 2018. 
0.177 
0,128 
0.095 
0.072 
0,055 
0,043 
0.035 
0.028 
0,023 
0,019 
0.016 
0,014 
0.012 
0.011 
0,009 


_NIKFAC 
12.862 
Ti 2T2 
“42,296 
7325269 
“9,362 
“6.553 
ee! 
70.946 
Otge 
0.468 
0.666 
0.859 
0,955 
0.986 
0.996 
0.998 
0.999 
1.000 
14.000 
1.000 
1,000 
1.000 
1.000 
12000 
1,000 
1,000 
peerayea 
1,000 


188 











Cc 


O\ZE;s; ZARB 8, 8B y SY yarawag) 


[xGeA © YTIDOISY Js AMASHTORT ow 








; 
for,0 ©&@ | $B\W OELaE -\e 8, TIGER 
‘N99 Cf SleaVVAOY ZADOW 2, TE=@ SYRRG 
me UAZUOSAY GAM £. ETLeQuyos QuHODTS 
“0 .@ 2i TIS AQG ZOLTALIVaES CRAGHATS 
t oWtwS,O,Utxd, 2 QHOTIAT NOLTAKAISA 
oe CUA 2¢ WREWRSE BI ROLTAT TMETAO 
WOSBAGATHISHOLTAAUSIAD SAO TERRE 
<a GSTTA Lo JVeLt JV qMat 
cr EE seeve..9 £.0 
ta.0 eeeee.0 
t TST .0 . teeee.o 
rr ete .o crece.o 
3 ' seri Teece .O 
ce 7 P&E f r¥veee.o 
taé.e€ pLe.r Gzeee ,O 
v6.0 rov.i selee.o 
er.t Ber: saeee .o 
’ ~ v8. eveeec.o 
0 edz. 0 eaeee .0 
Te eee ,0 
sé 4 rés.. © weeece.o 
Be .6 TTs 0 ##eee.0 
) Bif.0 réeee.o 
6ee,0 €e&0,0 *iece,.o 
0 St@6.0 2e8ee.0 
7 cco ,0 sTeee.o 
6¥9.0 pveee.o 
¢ eréee.o 
csvee.o 
: eevte,o 
eeseec .o 
evdee.0 
#e2zee.0 
steer, 0 
p#E Oo. 








a) 





3° 


- 


EXPERIMENTAL DATA FOR CRYSTAL I2 
VCTIJ+VvC0J 


TEMP 
2.966 
2.955 
2.930 
2,897 
2.867 
2.811 
2.766 
22724 
2.693 
2.664 
2.635 
2.604 
2.580 
2;552 
24522 
2.480 
2.434 
2,366 
2,279 
2.237 
2.189 
2.147 
2.103 
061 
1.983 
1.934 
1.873 
1.988 
1.935 
1.896 
1.867 
1.756 
1.665 
1.596 
1.525 
1.449 
e373 
1.264 
1207 
1.088 
1.008 
0.926 
0.648 


vir 
663.50 
663,55 
663,64 
663,80 
663,94 
664,18 
664,39 
664,56 
664,70 


664,78 © 


664,89 
604,99 
605.08 
665.17 
665, 26 
665,33 
b65.52 


BOOS 67 


605.91 
666,00 
660,11 
666.20 
666,28 
666.33 
666,46 
600,54 
666.61 
666,70 
666.76 
666,79 
666, 81 
666,90 
666,97 
666,993 
667.903 
667.06 
667.093 
OCW? 
667,13 
667,16 
667,16 
667,16 
607,16 


0,99421 
0.99429 
0.99442 
0.39466 
0.939487 
O,ago2s 
0.99554 
0.99580 
0.99601 
0.99613 
0.99629 
0.99644 
0.99658 
O7996 70 
0.99684 
0.93703 
0.99723 
0.99754 
0.99782 
0.99796 
0.99812 
0.99826 
0.99837 
0.99845 
0.99864 
0.33876 
0.99886 
0.99907 
0.99909 
0.99913 
Os393 17 
0.99930 
0.993940 
0.99944 
0.99950 
0.99954 
0.93959 
0.99963 
0.99965 
0.999693 
0.399969 
0797209 
0.99969 


NIMPAC 
1.008 
1.008 
1.002 
4004 
1.001 
0,998 
1.000 
1,004 
1.010 
0.998 
0,998 
0.995 
0.998 
0.994 
0.992 
0.991 
0.988 
0.985 
0,981 
0.976 
0.975 
0.975 
0,967 
0.951 
0.937 
0,937 
0,921 
45125 
1.116 
15100 
1.095 
1,056 
1.032 
0.977 
0.941 
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